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The  input  impedance  of  a  transistor  is  usually  difficult  to  predict  either 
by  compact  models  or  two-dimensional  (2D)  device  simulations  due  to  the 
three-dimensional  (3D)  distributed  effect  in  the  input  port  of  the  device.  The 
voltage  and  current  distribution  characteristics  vary  at  different  biases  as  well 
as  at  different  operating  frequencies  due  to  this  effect  resulting  in  the  modula- 
tion of  input  characteristics.  In  this  dissertation,  an  analytical  model 
describing  the  DC  and  AC  voltage  and  current  distributed  effects  in  the  polysil- 
icon  and  intrinsic  base  regions  of  stripe  bipolar  junction  transistors  (BJTs) 
with  double  polysilicon  technology  is  presented.  It  is  shown  that  the  bias  and 
frequency  dependent  distributed  effect  in  the  base  polysilicon  contacts  cause 
an  unequal  division  of  base  current  between  the  two  base  polysilicon  contacts 
and  result  in  a  re-distribution  of  base  current  in  the  base  regions  at  different 
levels  of  current  injection  and  operating  frequencies.  Therefore,  the  input 
characteristics  are  also  modulated  at  different  operating  conditions.  It  is  also 
observed  that  this  effect  has  different  degrees  of  influence  on  devices  with 


VI 


various  types  of  layout  structures,  geometry  sizes,  and  polysilicon  contact  resis- 
tances. Thus,  resulting  performance  of  these  devices  is  quite  different.  The 
distributed  effect  on  the  transient  and  AC  performance  of  advanced  metal- 
oxide-semiconductor  field-effect  transistors  (MOSFETs)  is  also  studied  based 
on  the  3D  device  simulations,  and  the  results  show  the  distributed  effect  is  less 
important  for  MOSFETs. 
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CHAPTER  1 
INTRODUCTION 


Modern  high  speed  BJTs  are  built  with  highly  scaled  lateral  and  vertical 
device  geometries.  The  vertical  scaling  increases  the  current  densities  and 
reduces  the  minority  carrier  storage  whereas  the  lateral  scaling  decreases  the 
parasitic  capacitances  and  resistances  in  the  devices. 

In  single  polysilicon  bipolar  transistors  (refer  to  Fig.  1.1(a)),  the  heavily 
doped  polysilicon  layer  is  used  as  a  diffusion  source  to  form  a  shallow  emitter 
junction  and  serves  as  a  low  resistivity  contact  for  the  emitter.  The  polysilicon 
emitter  structure  has  also  allowed  scaling  of  the  base  width  while  simulta- 
neously improving  the  device  current  gain  and  performance  [War95].  The 
double-polysilicon  self-aligned  (DPSA)  bipolar  device  structure  (refer  to  Fig. 
1.1(b))  first  appeared  over  ten  years  ago  [Nak80,  Tan82]  and  became  the  key 
technology  for  most  high  performance  BJTs  [Yam93,  Ker92,  Tor91].  The  use  of 
a  heavily  doped  polysilicon  layer  as  the  base  contact  scales  the  extrinsic  base 
regions  and  results  in  much  smaller  parasitic  capacitance  compared  with  the 
single  polysilicon  transistors  [Nak95].  The  double  polysilicon  processes  with 
sidewall  spacer  technologies  reduce  the  link  between  the  intrinsic  and  extrinsic 
base  regions  and  decrease  the  base  resistance  [Zde87]. 

The  resistivity  in  the  base  region  causes  a  debiasing  effect  along  the 
distributed  base  current  path  and  results  in  a  nonuniform  voltage  and  thus 
emitter  current  injection  across  the  base-emitter  junction.  The  emitter  current 
density  at  the  edge  is  higher  than  that  at  the  center  of  the  intrinsic  base  region 
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Figure  1.1  The  cross  section  of  (a)  a  single  polysilicon  and  (b)  a  double 
polysilicon  BJTs.  Note  the  longer  inter-connection  between 
intrinsic  and  extrinsic  base  regions  in  the  single  polysilicon  BJT. 


which  is  the  so  called  "current  crowding"  effect.  Current  high  speed  BJTs  tend 
to  operate  at  high  current  density  and  high  frequency  where  the  current  crowd- 
ing effect  is  prominent.  The  effect  in  the  intrinsic  base  region  on  the  DC  and 
AC  base  resistance  has  been  widely  investigated  [Yua88,  Rhe95,  Gho65,  Ver91, 
Cho90].  The  lateral  scaling  of  the  emitter  stripe  width  with  DPS  A  technologies 
alleviates  the  effect.  The  base  polysilicon  contact  is  also  important  for  its  role 
in  determining  the  base  resistance  and  a  similar  debiasing  effect  takes  place  in 
the  polysilicon  region  as  well.  As  device  dimensions  are  scaled,  the  resistance 
through  the  base  polysilicon  contact  increases  since  the  area  of  the  polysilicon 
layer  must  also  be  scaled.  Thus,  the  debiasing  effect  in  the  polysilicon  base 
region  is  even  more  important  in  the  scaled  DPSA  transistors. 

Contemporary  MOSFET  technologies  utilize  polysilicon  as  the  gate 
material  to  reduce  gate  (channel)  length  as  well  as  control  threshold  voltage. 
Usually  the  metal  contacts  of  the  polysilicon  gate  are  attached  at  the  end  of  the 
gate.  The  distributed  voltage  and  current  effects  similar  to  that  of  BJTs  are 
observed  in  the  gate  of  MOSFETs  at  high  speed  operations.  Shrinking  the  gate 
length  greatly  improves  the  speed  and  increases  current  density.  The  thickness 
of  the  gate  is  correspondingly  scaled  due  to  the  required  planarity  Both  of  the 
tendencies  toward  high  speed  operation  and  the  scaling  of  feature  sizes  in 
MOSFETs  enhance  the  distributed  characteristics  of  the  device. 

Two-dimensional  device  simulations  and  compact  model  simulations 
currently  are  the  main  tools  to  characterize  the  performance  of  the  semicon- 
ductor devices.  However,  the  nonuniformity  of  voltage  and  current  in  the 
polysilicon  contact  region  cannot  be  properly  described  by  these  methods  due 
to  the  distributed  voltage  and  current  characteristics  that  occur  in  polysilicon. 
Therefore,  the  understanding  and  modeling  of  the  distributed  effects  is  impor- 
tant for  accurate  and  efficient  design  of  high  speed  semiconductor  devices.  The 


goals  of  this  work  are  to  (i)  characterize  the  distributed  and  debiasing  effects, 
(ii)  model  the  effects,  and  (iii)  evaluate  the  impact  of  the  effects  on  the  DC  and 
AC  device  performance. 

The  device  simulator  FLOODS  (FLorida  Object-Oriented  Device 
Simulator)  [Lia94a,  Lia94b,  Lia94c]  is  the  main  thrust  for  this  study.  To 
accurately  describe  the  complicated  physical  mechanisms  of  modern  semicon- 
ductor devices,  several  new  models  and  functions  have  been  added  to  FLOODS. 
The  development  and  implementation  of  these  functions  are  presented  in 
chapter  2. 

In  chapter  3,  a  quasi-3D  DC  distributed  model  which  considers  the  effect 
in  both  the  intrinsic  and  polysilicon  base  regions  is  developed  for  BJTs  and 
applied  to  the  computation  of  base  resistance.  Model  verification  is  accom- 
plished through  comparison  with  the  measured  base  resistance. 

In  chapter  4,  a  quasi-3D  AC  distributed  model  is  developed  to  analyze 
the  current  redistribution  in  BJTs  at  different  operating  frequencies.  The 
model  is  also  applied  to  FLOODS  for  quasi-3D  simulations  in  which  nonunifor- 
mity  of  voltage  and  current  in  the  third  dimension  are  taken  into  account.  The 
model  is  verified  by  the  comparison  between  the  simulations  and  s-parameter 
measurements. 

In  chapter  5,  a  3D  device  simulation  scheme  and  device  structures  are 
presented,  and  the  problems  associated  with  the  3D  simulations  are  discussed. 
The  DC  distributed  effects  presented  in  BJTs  are  studied  with  3D  simulations 
and  compared  with  the  results  of  the  3D  model.  The  layout  structure,  geometry 
size  and  polysilicon  sheet  resistance  dependence  of  the  BJT  performance  due 
to  the  distributed  effects  are  also  discussed. 

In  chapter  6,  AC  distributed  effects  in  BJTs  are  studies  using  3D  simula- 
tions. The  simulation  results  confirm  the  accuracy  of  the  AC  model  derived  in 


chapter  4.  The  AC  performance  of  devices  with  different  layout  structures, 
geometry  sizes  and  polysilicon  sheet  resistances  are  then  characterized  at 
different  frequencies.  The  impacts  of  the  distributed  effects  on  the  device  and 
circuit  design  of  high  speed  B  JTs  are  also  presented. 

The  implementation  of  a  transient  simulation  scheme  in  FLOODS  is 
presented  in  chapter  7.  The  effects  of  the  distributed  gate  resistance  on  the 
transient  and  AC  characteristics  of  a  0.18um  NMOS  transistor  with  different 
gate  sheet  resistances  are  studied  based  on  the  results  of  3D  device 
simulations. 


CHAPTER  2 

IMPLEMENTATION  OF  NEW  PHYSICAL  MODELS 

AND  FUNCTIONS  IN  FLOODS 


2.1  Introduction 

FLOODS  (FLorida  Object  Oriented  Device  Simulator)  [Lia94b,  Lia94c] 
is  a  general  purpose  3D  device  simulator.  In  addition  to  Poisson  and  two  carrier 
continuity  equations,  lattice  self-heating  and  two  carrier  energy  balance 
equations  can  be  solved  in  FLOODS  to  predict  the  characteristics  of  advanced 
devices  with  self-heating  and  hot  carrier  effects.  Unlike  MEDICI  which  de- 
couples the  energy  balance  equations  from  the  carrier  continuity  equations, 
FLOODS  solves  the  six  equations  simultaneously  to  achieve  self-consistency  in 
these  equations.  Another  advantage  of  FLOODS  is  that  it  is  implemented  in 
C++  which  provides  benefits  of  modularity,  great  portability,  and  reusability  of 
code  modules. 

C++  and  object-oriented  programming  techniques  are  characterized  by 
the  inheritance  and  dynamic  binding  functions.  Inheritance  allows  a  new  class 
to  inherit  the  basic  functions  and  data  of  its  parents.  No  additional  program- 
ming is  required  in  the  new  class  except  for  those  functions  or  data  that  either 
extend  or  replace  characteristics  inherited  from  their  parents.  Dynamic 
binding  is  provided  by  the  virtual  functions:  the  program  resolves  which 
virtual  function  in  the  inheritance  hierarchy  to  call  based  on  the  run-time  type 
of  the  actual  object.  Run- time  type  resolution  encapsulates  the  implementation 
details  of  an  inheritance  hierarchy  from  the  programs  that  make  use  of  it. 


Users  do  not  have  to  know  how  the  objects  are  constructed  and  the  user  code  is 
not  subject  to  change  to  the  hierarchy.  These  two  features  greatly  simplify  the 
implementation  of  new  models  or  functions  in  FLOODS. 

In  this  chapter,  a  single-node  self-heating  model  is  developed  in 
FLOODS  to  reduce  the  running  time  of  full  thermodynamic  simulations 
without  losing  much  accuracy.  The  implementation  of  the  current  boundary 
condition  is  shown  which  provides  unique  and  stable  solution  at  the  situations 
of  small  differential  resistance  or  negative  resistance.  A  carrier  tunneling 
model  in  the  polysilicon  emitter  interface  is  exhibited  which  gives  a  proper 
description  of  the  carrier  distribution  and  conduction  in  the  boundary.  The 
coupling  of  a  circuit  simulator  with  FLOODS  is  also  presented  to  include  the 
parasitic  effects  in  device  simulations  and  provide  direct  circuit  simulations 
with  complex  device  structures. 

2.2  Single  Node  Self-Heating  Model 
2.2.1  Motivation 

The  traditional  drift- diffusion  (DD)  transport  model  in  device  simula- 
tions solves  three  semiconductor  equations:  Poisson  equation  plus  the  electron 
and  hole  continuity  equations.  However,  as  device  technologies  improve,  these 
three  equations  become  insufficient  to  describe  the  characteristics  of  state-of- 
the-art  semiconductor  devices.  The  failure  comes  from  the  fact  that  a  quasi- 
static  condition  is  assumed  while  deriving  the  DD  transport  equations  which  is 
not  valid  for  carrier  transport  with  rapid  spatial  and  time  variation  of  the  local 
electric  field.  The  large  field  and  field  gradient  inside  the  scaled  devices  cause 
lattice  self-heating  and  hot  carrier  nonlocal  transport  effects.  Lattice  self- 
heating  and  carrier  energy  balance  equations  [Str62,  Coo83,  Che93,  Fus95b] 
govern  the  energy  transfer  among  electrons,  holes,  and  the  crystal  lattice  which 
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take  into  account  most  of  the  possible  mechanisms  associated  with  nonlocal 
and  self-heating  effects.  The  full  thermodynamic  (TD)  simulations  solve  these 
three  energy  equations  together  with  DD  equations  producing  better  predic- 
tion capabilities  compared  with  the  DD  simulations  [Lia94b,  Ago94] . 

The  specific  position  of  a  simulated  device  geometry  is  represented  by  a 
node  in  the  mesh.  Equations  are  solved  and  solutions  are  stored  at  each  node. 
For  a  simulated  structure  of  N  nodes,  the  total  equations  M  solved  in  a  TD 
simulation  are  6xN  whereas  3xJV  equations  are  solved  in  a  DD  simulation. 
The  required  memory  and  CPU  time  to  solve  the  problem  are  approximately 
proportional  to  M1,5.  Thus,  the  memory  and  CPU  time  requirements  of  TD 
simulations  are  2.8  times  of  the  DD  simulations  for  the  same  number  of  nodes. 
A  scaled  device,  in  general,  requires  thousands  of  nodes  to  define  the  device 
structure.  Although  TD  simulations  are  more  accurate,  the  simultaneous 
solving  of  six  equations  is  CPU  and  memory  intensive  compared  with  DD 
simulations. 

In  this  subsection,  an  algorithm  which  reduces  the  solution  time  of  the 
TD  simulations  with  little  sacrifice  in  the  accuracy  is  presented. 

2.2.2  Model  Development 

The  energy  balance  equations  for  electrons  and  holes  are  given  by 

li(n.Un)  =  -VSn-Rwn  (2-1) 


and 


§-t(p-Up)  =-V$p-Rwp  (2-2) 


where  n  and  p  are  the  electron  and  hole  densities,  respectively,  Un  „  is  the 

wp 


average  local  carrier  energy,  Sn,  p  is  the  carrier  energy  flux,  and  Rwn  and  R 
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are  the  energy  transfer  rates  from  electron  and  hole,  respectively,  to  other 
subsystems  due  to  the  scattering  or  carrier  recombination/generation.  The 
lattice  self-heating  equation  governing  joule-heating  and  heat  conduction  in 
the  semiconductor  material  is 

^(Ul)=-V*Sl-Rwl  (2-3) 

where  U^  is  the  local  lattice  energy,  Sl  is  the  lattice  heat  flux,  and  Rwl  *s  ^ne 
local  energy  transfer  rate  from  crystal  lattice  to  the  carrier  subsystems.  These 
three  equations  basically  state  the  conservation  of  energy  at  a  specific  position 
of  the  device  (node)  during  a  period  of  time.  If  no  external  energy,  for  example, 
photon,  is  supplied,  the  complete  thermodynamic  system  is  energy  conserva- 
tive and  we  have  RWL  +  Rwn  +  R  =  0  .  The  energy  conservation  relation  and 
the  energy  flow  between  the  three  subsystems  are  plotted  schematically  in 
figure  2.1.  The  discretization  scheme  of  eq.  (2-3)  with  box  (finite  difference) 
method  is  shown  in  figure  2.2  and  can  be  expressed  by 

l(UL)'Ai  =  -J,SLb-lb-^L-lj-RWL-Ai  (2-4) 

b  j 

where  A  j  is  the  box  area  of  a  node  i,  L  is  the  length  of  the  box  edge  between  node 

i  and  bulk  node,/,  and  1%  is  the  length  of  box  edge  between  node  i  and  boundary 

b.  The  boundary  condition  of  lattice  self-heating  is  given  by 

2^SLb-lb  =  Gth-(TLi-TL{amhi(int^)  (2-5) 

b 

where  Gt^  is  the  thermal  conductance  of  the  interconnection  boundary  and  T^i 

is  the  lattice  temperature  at  node  i. 

Sl  in  the  bulk  semiconductor  is  primarily  the  heat  flow  due  to  the  lattice 

temperature  gradient  and  can  be  expressed  as 
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f 

electron  subsystem 


Figure  2.1.  Schematic  diagram  of  the  thermodynamic  system.  The  energy 
system  are  divided  into  three  subsystems:  electrons,  holes  and 
lattice.  Energy  in  each  subsystem  can  transfer  to  the  other  two 
subsystems  by  way  of  scattering,  recombination  or  generation. 
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Figure  2.2.  The  discretization  scheme  with  a  box  method  in  a  two  dimensional 
mesh.  The  shaded  area  indicates  the  area  A^  associated  with  node 
i.  In  single  node  self-heating  model,  T^  ;  =  TV/- 
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Sl  =  -KL  •  VTL  (2-6) 

where  /fL  is  the  lattice's  thermal  conductivity.  The  active  region  of  a  scaled  BJT 
or  MOS  transistor  is  usually  small  and  the  local  self-heating  of  the  silicon 
lattice  can  be  distributed  over  the  active  region.  If  we  assume  the  lattice 
temperature  is  constant  in  the  bulk  region,  we  have  VTL  =  0  and  Sl  =  0  and 
eq.  (2-4)  becomes 

A([/L) .  ^  =  -JSLb  ■  lb-RWL  ■  At  (2-7) 

b 

Since  the  lattice  temperature  is  constant  in  the  bulk,  the  solution  of  T^  in  one 
node  can  represent  T/^s  at  the  other  N-l  nodes  and  only  one  instead  of  N  lattice 
self-heating  equations  needs  to  be  solved  in  the  bulk  mesh  of  N  nodes.  The  size 
of  the  Jacobian  matrix  becomes  (5iV +  1)X  (5N  +  1)  whereas  the  TD  simula- 
tion is  required  to  solve  a  6N  x  6N  matrix.  Though  solving  one  equation, 
FLOODS  loops  over  the  N  nodes  and  sums  the  energy  transfer  rate  Rwl  of  each 
node  in  the  same  equation.  Thus,  eq.  (2-7)  can  be  interpreted  as  follows:  the 
total  lattice  energy  generated  in  the  device  is  represented  by  a  single  temper- 
ature in  one  node  and  balanced  by  the  heat  flux  at  the  contact  boundaries.  In 
the  isothermal  model,  the  self-heating  equation  is  not  taken  into  account  and 
the  size  of  the  Jacobian  matrix  is  5N  x  5N  which  is  about  the  same  size  as  the 
matrix  in  the  single  node  self-heating  (SNSH)  model  if  N  »  1. 

2.2.3  Results  and  Discussion 

The  lattice  temperature  distribution  in  an  npn  BJT  with  TD  simulation 
of  FLOODS  is  shown  in  figure  2.3.  The  base-emitter  junction  is  forward  biased 
at  1.0V  which  drives  the  device  into  the  high  current  injection  regime  and 
pushs  the  base-collector  space-charge  region  to  the  edge  of  the  collector  buried 
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Figure  2.3.  The  lattice  temperature  contour  in  an  npn  BJT  at  the  biases  of  Vbe 
=  1.0V  and  Vbc  =  3V.  The  ambient  temperature  is  300K  and  the 
thermal  conductances  at  the  contacts  are  infinite.  Dark-shaded 
area  is  the  emitter  and  light-shaded  area  is  the  base. 
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layer  where  the  hot-carrier  effect  takes  place  due  to  the  high  electric  field.  The 
electron  temperature  in  the  region  is  as  high  as  4500K  and  the  energy  transfer 
from  electron  to  other  subsystems  is  significant.  However,  the  peak  lattice 
temperature  in  the  region  is  only  6K  higher  than  the  boundaries.  The  base  and 
collector  currents  are  mainly  determined  by  the  drift-diffusion  in  the  emitter 
and  base  regions,  respectively,  and  the  self-heating  effect  shows  up  in  the 
physical  models  such  as  intrinsic  carrier  density  and  mobility.  Since  the 
temperature  distribution  is  uniform  in  these  regions,  the  constant  temperature 
assumption  is  valid  in  the  B  JT  devices.  The  comparisons  of  the  running  time 
and  error  between  these  three  simulations  are  summarized  in  Table  2.1.  At  the 
biases  of  Vbe=0.9V  and  Vbc=3V,  the  error  of  the  SNSH  model  is  less  than  2% 
but  the  solution  time  is  20%  less  than  the  TD  simulation,  whereas  the  error  of 
the  isothermal  simulation  is  greater  than  7%  and  running  time  is  slightly 
faster  than  the  SNSH  simulation. 

Table  2.1  The  comparisons  among  three  simulation  schemes  for  BJTs 


Model 

Vemit  =  -0.8V,  Vcoll  =  3.0V 

Vemit  =  -0.9V,  Vcoll  =  3.0V 

CPU  time 

error  (%) 

CPU  time 

error  (%) 

Full-thermo- 
dynamic 

1 

0 

1 

0 

Single  Node 
Self-Heating 

0.82 

0.10 

0.77 

1.52 

Isothermal 

0.78 

1.17 

0.71 

7.09 

The  lattice  temperature  along  the  surface  of  a  0.6jim  NMOS  with  TD 
simulation  is  plotted  in  figure  2.4  when  the  biases  at  drain  and  gate  are  7V.  Due 
to  the  thermal  conductances  at  all  contacts,  the  device  is  heated  up  and  the 
lattice  temperature  at  the  boundaries  is  at  least  40K  higher  than  the  ambient 
temperature.  At  the  same  operation  condition,  the  lattice  temperature  of 
SNSH  and  isothermal  simulations  is  345K  and  300K,  respectively  Since  the 
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Figure  2.4.  Lattice  temperature  along  the  surface  of  a  NMOS  from  TD  and 
SNSH  simulations.  Both  the  biases  at  gate  and  drain  are  7V  and 
the  ambient  temperature  is  300K.  A  thermal  conductance  of  9000 
K/Wcm  is  attached  to  each  contact.  The  arrows  mark  the  gate 
width. 
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drain  current  is  governed  by  the  electron  transport  in  the  channel  and  the  self- 
heating  effect  shows  up  in  the  mobility  model,  the  large  difference  of  channel 
lattice  temperature  between  isothermal  and  TD  simulations  results  in  signifi- 
cant error  in  drain  current.  In  contrast,  the  discrepancy  of  drain  current 
between  SNSH  and  TD  simulations  is  little  due  to  the  small  difference  of  the 
channel  lattice  temperature.  Figure  2.5  shows  the  comparisons  among  the  TD, 
SNSH,  and  isothermal  simulations  of  the  drain  current  as  a  function  of  biases. 
The  comparison  between  TD  simulations  and  experimental  data  is  shown  in 
Liang  and  Law  [Lia94b].  The  drain  current  predicted  by  the  TD  model  is 
always  lower  than  that  of  the  isothermal  model  because  the  channel  lattice 
temperature  in  TD  simulation  is  higher  than  that  of  the  isothermal  simulation 
which  increases  the  lattice  scattering  and  decreases  the  electron  mobility.  The 
degradation  of  the  electron  mobility  predicted  by  TD  model  causes  negative 
output  conductance  at  high  Vgs  which  is  not  observed  in  the  isothermal  simula- 
tion. The  error  of  the  isothermal  simulation  is  larger  than  10%  when  Vgs  =  Vds 
=  7V.  For  the  SNSH  simulation,  the  degradation  of  mobility  and  the  negative 
conductance  are  also  observed.  The  error  of  SNSH  simulation  is  within  1.2% 
but  the  simulation  time  is  20%  less  than  the  TD  simulation.  The  cpu  time  and 
error  of  the  three  simulations  are  summarized  in  Table  2.2. 

The  accuracy  of  this  approach  is  determined  by  the  difference  between 
the  lattice  temperature  of  SNSH  simulation  and  actual  lattice  temperature  in 
the  critical  active  regions,  i.e.  channel  in  MOSFETs  and  intrinsic  base,  emitter 
and  base-collector  space  charge  regions  in  BJTs.  For  a  device  with  large 
temperature  variations  in  a  big  critically  active  area  such  as  a  power  transistor, 
the  physical  mechanisms  can  not  correctly  be  described  by  a  single  lattice 
temperature.  This  could  result  in  a  larger  error  in  the  terminal  current. 
However,  this  approach  appears  useful  for  a  submicron  device  structure  where 
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the  lattice  temperature  variation  in  the  critically  active  region  is  small  and  can 

be  properly  approximated  by  the  SNSH  simulations.  The  accuracy  of  the  SNSH 

model  is  insensitive  to  the  size  of  the  noncritically  active  region  since  there  is 

no  heat  generated  in  the  area  and  it  has  little  effect  on  the  lattice  temperature 

of  the  SNSH  simulations. 
Table  2.2  The  comparisons  among  three  simulation  schemes  for  MOSFETs 


Model 

Vdrain  = 

7V,  Vgate  =  6V 

Vdrain  = 

7V,  Vgate  =  7V 

CPU  time 

error  (%) 

CPU  time 

error  (%) 

Full-thermo- 
dynamic 

1 

0 

1 

0 

Single  Node 
Self-Heating 

0.80 

1.02 

0.72 

1.16 

Isothermal 

0.76 

8.69 

0.78 

10.5 

2.3  Current  Boundary  Condition 

In  the  analysis  of  BJT  breakdown  or  CMOS  latch-up  triggering,  an  I-V 
curve  snake  back  phenomenon  is  usually  observed.  The  current  is  a  multi- 
valued function  of  the  applied  voltage  in  that  regime  which  implies  that  the 
numerical  solutions  of  the  analysis  may  end  up  with  a  solution  in  one  of  three 
distinct  and  stable  states  depending  on  the  initial  guess.  Furthermore,  in  the 
negative  resistance  regime  or  around  the  snake  back  point  where  dV/dl  =  0,  it 
is  very  difficult  or  impossible  to  obtain  a  solution  with  a  voltage  boundary 
condition.  Since  the  current  is  a  single-valued  function  of  the  applied  voltage 
in  these  regimes,  the  current  boundary  condition  can  provide  a  unique  and 
faster  converged  solution  than  the  voltage  boundary  condition. 

For  device  simulators  such  as  MEDICI,  the  boundary  condition  for  the 
Poisson  equation  in  the  contact  region  is  to  fix  the  potential.  The  definition  is 
invalid  if  there  is  carrier  accumulation  or  depletion  at  the  contact  where  the 
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relative  position  of  the  band-gap  and  majority  carrier  Fermi  energy  begins  to 
change.  In  FLOODS,  the  voltage  boundary  condition  for  Poisson  equation  is 
defined  by  the  fixed  Fermi-level  EFC  at  the  contact  and  given  by  [Lia94a] 

EFC{y,  C,  Tc)  =  EFC0(yo,  C0,  Tco)  -  Vapplied  (2-8) 

where  V|/,  C  and  Tc  represent  the  potential,  majority  carrier,  and  majority 
carrier  temperature  solutions,  respectively.  Subscript  o  indicates  the  equilib- 
rium solution  variables  and  VappUed  is  the  applied  voltage  at  the  contact.  For 
the  current  boundary  condition,  eq.  (2-8)  is  replaced  by 

JappUed  =  X(.Jni^n,Tn)+Jpi{Xf,p,Tp)+JdiSii(^n,p)) 

i 

(2-9) 

where  Jappued  is  tne  applied  current  density  at  the  contact,  i  indicates  the 
contact  nodes,  n  andp  are  electron  and  hole  solutions,  respectively,  and  Jni>  Jpi, 
and  Jji.i  are  the  electron,  hole,  and  displacement  current  densities, 
respectively. 

The  dependence  of  current  densities  on  potential  \\f  is  weaker  than  on 
Fermi  energy.  Thus,  in  the  nonlinear  Newton  iteration,  the  convergence  of  the 
current  boundary  condition  is  slower  than  that  of  the  voltage  boundary  if  the 
variation  of  \|/  solution  is  large,  i.e.  dl/dV  is  small.  However,  the  dependence  of 
the  current  densities  on  carrier  density  is  slightly  stronger  than  the  Fermi 
energy.  Therefore,  the  convergence  of  the  current  boundary  will  be  faster  than 
that  of  the  voltage  boundary  if  the  current  \\f  solution  is  very  close  to  the  previ- 
ous one  and  the  variation  of  the  carrier  solution  is  large  in  the  Newton 
iteration,  i.e.  dl/dV  is  large.  The  Newton  iteration  steps  required  for  two- 
dimensional  BJT  simulations  with  voltage  and  current  boundary  conditions 
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Figure  2.6.  The  number  of  Newton  iteration  required  for  normal  active  BJT 
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are  shown  in  figure  2.6  as  a  function  of  collector  current  density.  It  can  be  seen 
that  the  voltage  drive  boundary  condition  is  preferred  at  low  current  regime 
where  dl/dV  is  small.  In  contrast,  the  current  drive  boundary  has  quicker 
convergence  in  the  high  current  regime  where  dl/dV  is  large. 

2.4  Interface  Transport  Model  of  Polvsilicon  Emitter  Contacts 
2.4.1  Motivation 

By  the  use  of  heavily  doped  polysilicon  as  the  diffusion  source  of  emitters 
and  emitter  contacts,  the  fabrication  of  self-aligned,  shallow  junction  bipolar 
transistors  has  been  made  possible  and  it  leads  to  the  features  of  higher 
current  gain  [Gra76,  Gra79],  packing  densities,  and  much  higher  frequency 
capabilities  [Suz85,  ParS6,  Ich89]. 

A  polysilicon  layer  is  composed  of  mono-crystalline  grains  and  disor- 
dered grain  boundaries.  Due  to  the  potential  barrier  and  recombination  in  the 
grain  boundaries,  polysilicon  has  a  smaller  carrier  mobility  and  a  shorter 
carrier  life-time  compared  with  mono-crystalline  silicon  [Set75].  A  reduced 
mobility  model  in  the  polysilicon  region  has  been  proposed  to  explain  the 
current  gain  enhancement  in  these  devices  based  on  the  retarded  transport 
mechanisms  of  minority  carriers  [Nin80].  During  the  fabrication  process,  a 
thin,  oxide-like  insulator  is  presented  at  the  mono-crystalline  silicon/polysili- 
con  interface  and  this  layer  has  a  higher  energy  band-gap  than  silicon  which 
acts  as  a  potential  barrier  to  carriers.  Graaff  [Gra79]  proposed  a  carrier  tunnel- 
ing model  to  explain  the  carrier  transport  through  the  oxide-like  layer  and  the 
current  gain  enhancement.  These  two  mechanisms  have  been  combined  in  the 
numerical  models  [Elt82,  Yu84]  which  describe  the  experimental  I-V  character- 
istics more  satisfactorily.  Chor  [Cho93]  further  applied  the  combination  model 
to  analyze  the  dependence  of  high  frequency  performance  on  the  polysilicon 
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Figure  2.7.  (a)  The  mesh  structure  in  a  contact.  The  boundary  conditions  are 
implemented  inside  the  interface  nodes  and  the  solution  variables 
of  each  node  are  the  flux  came  from  the  corresponding 
semiconductor  node,  (b)  The  band  diagram  of  a  polyemitter  contact. 


23 

emitter  transistors  and  showed  a  significant  difference  in  the  charge  storage 
behaviors  with  and  without  the  presence  of  polyemitter  contact. 

Man}^  parameters  associated  with  the  properties  of  the  grain  and  grain 
boundary  in  the  polysilicon  material  are  required  in  the  retarded  mobility 
models.  However,  the  physical  properties  of  the  polysilicon  are  not  well  under- 
stood and  vary  with  different  process  technologies.  Fossum  [Fos80]  proposed  to 
characterize  the  minority-carrier  behavior  in  polysilicon  with  an  effective 
mobility  and  a  lifetime  constant,  thus  resulting  in  a  reduced  effective  surface 
recombination  velocity  (ERV)  as  viewed  from  the  monocrystalline  emitter.  The 
concept  of  ERV  greatly  simplifies  the  boundary  condition  used  in  the  polyemit- 
ter contacts  of  BJTs  in  the  device  simulations.  Though  the  use  of  constant  ERV 
well  characterizes  the  base  current  at  low  and  intermediate  bias,  the  current 
at  high  biases  is  not  properly  described  due  to  the  nonlinear  I-V  characteristics 
of  the  interface  tunneling  current  [YU84] .  Thus,  the  boundary  condition  with 
the  combination  of  ERV  and  the  current  tunneling  model  is  used  in  FLOODS 
to  appropriately  represent  the  properties  of  a  polyemitter  contact. 

2.4.2  Model  Implementation 

To  simplify  the  implementation,  no  additional  material  or  mesh  is 
created  to  adopt  the  current  tunneling  model.  The  model  is  directly  applied  to 
the  boundary  conditions  with  the  current  algorithm  of  FLOODS.  Figure  2.7(a) 
illustrates  the  boundary  contact  structure  and  the  solution  variables  associ- 
ated with  each  mesh  in  FLOODS.  The  governing  equations  of  the  boundary 
conditions  are  solved  in  the  interface  mesh  where  the  solution  variables  are  the 
flux  O  of  potential  and  carriers.  Because  the  fixed  potential  boundary  condition 
is  invalid  in  the  presence  of  charge  accumulation  or  depletion  in  the  contacts, 
FLOODS  uses  a  fixed  Fermi  energy  boundary  condition  and  assumes  the  varia- 
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tion  of  the  majority  carrier  Fermi  energy  EFis  corresponding  to  the  change  of 
potential  with  opposite  sign,  i.e.,  dEF/dy  =  -l.  Thus,  the  boundary  condition  of 
Poisson  equation  at  the  polyemitter  interface  can  be  written  as 

%  +  ^f~  ■  (EF  +  VappUed)  =  0  (2-10) 

H  '  Lox 

where  q  is  the  electronic  charge,  eox  is  the  permittivity  and  tox  is  the  thickness 
of  the  oxide-like  insulator,  and  Ov  is  the  potential  flux  through  the  interface 
due  to  the  charge  accumulation  or  depletion  at  the  boundaries  of  the  interface. 
In  an  n-type  polysilicon  contact,  the  tunneling  current  flux  through  the 
insulator  is  given  by  [Yu84] 

AeK  *dnkTn  _dnqAV 

Jtn  = 5  "    •    }    j   um  ;  •  (1  -e  )  (2-11) 

q(dnk)2    sin(7cdBAr„) 


and 


-K         (  -&\ 

J 


kT„  0~bp 


*P  "  >^2%mp    \-dpkTp 


P2~P\-e 


V 


(2-12) 


where  p2  is  the  hole  concentration  at  the  monocrystalline  side  of  the  interface 
and  pj  at  the  polysilicon  side,  AV  is  the  voltage  across  the  interface  which 
equals  EFn2  +  VappUed  (refer  to  figure  2.7(b)),  k  is  the  Boltzmann  constant,  mm 
and  mp  are  the  effective  masses  of  electron  and  hole,  respectively,  and  A  is  the 
Richardson  constant  of  electrons 

4Km„qk 

A  -  f~  (2-13) 

h 

If  AV  is  less  than  the  barrier  height  of  the  interface,  the  coefficients  bc  and  dc 
are  given  by 

4nt 


bc  =  -jf-  fi™n.pVKp  (2-14) 
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Figure  2.8.  The  current  components  in  each  region  of  the  polysilicon  emitter 
contact  for  (a)  electron  and  (b)  hole.  Si  is  the  surface  recombination 
velocity  of  the  interface,  and  n0  and  p0  are  the  electron  and  hole 
densities  at  thermal  equilibrium,  respectively. 
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and 

2TCtnr      \2mn  ' 

d'-^hr  <2-15) 

where  h  is  the  Planck  constant,  c  represents  electron  n  or  holep,  and  V^  is  the 
average  energy  barrier  height  of  the  interface  for  electrons  or  holes. 

The  carrier  flux  components  in  each  region  of  the  polysilicon  emitter 
contact  are  plotted  in  figure  2.8.  The  excess  carriers  AC  1  in  the  polysilicon  side 
of  the  interface  can  be  related  to  the  carrier  flux  Qqi  °f  the  same  region  by 
ACj  =  Cj  -  C0  =  c  .  Then,  the  boundary  condition  for  electron  or  holes  in 
the  interface  mesh  is  given  by 

1±— iS-W.-.T    =  0  (2-16) 


ERV, 


and  the  continuity  equations  for  those  nodes  in  the  semiconductor  mesh  at  the 
contact  boundary  can  be  expressed  by 

^tc±Sic(C2^C02))  +  YjJc2  =  0  (2-17) 

where  "+"  for  holes  and  "-"  for  electrons  if  Jtc  >  0,  and  vice  versa. 

2.4.3  Results  and  Discussion 

The  comparison  of  the  simulated  Gummel  plot  of  a  polyemitter  BJT  with 
different  thickness  of  interface  layers  is  shown  in  figure  2.9.  Without  an  inter- 
face layer,  the  base  current  remains  a  straight  line  at  high  biases  and  the 
collector  current  base  push-out  effect  takes  place  at  V^  =  0.96V.  With  the  intro- 
duction of  the  interface  layer,  both  the  base  and  collector  current  deviate  from 
the  curves  without  the  presence  of  the  interface  layer  (ideal  curves).  The 
thicker  the  layer,  the  larger  the  current  deviation. 
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Figure  2.9.  The  simulated  Gummel  plot  with  different  interface  thicknesses. 
The  voltage  deviation  AV  =  RbIb  +  ^E^-E- 
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Figure  2.10.  The  hole  densities  in  the  emitter  of  a  BJT  with  different 
thicknesses  of  interface  layers.  The  BJT  is  biased  at  V^g  =  1.0V 
and  Vcb  =  2.0V.  The  interface  layer  locates  at  x  =  0.0|xm  and  base- 
emitter  junction  is  at  x  =  0.85|J.m. 
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There  are  two  major  effects  associated  with  the  presence  of  the  insulator 
interface  in  the  poly  emitter  contact.  The  first  one  is  that  the  finite  tunneling 
probability  limits  the  carrier  transport  through  the  interface  and  results  in  an 
accumulation  of  holes  and  depletion  of  electrons  at  the  mono-crystalline 
emitter  side  of  the  interface  in  an  npn  BJT  A  thicker  interface  causes  more 
hole  accumulation  and  electron  depletion  due  to  the  smaller  tunneling  proba- 
bility. However,  the  collector  current  is  governed  by  the  characteristics  of 
electrons  in  the  base  region  and  is  not  affected  by  the  electron  depletion  in  the 
emitter  interface  region.  On  the  other  hand,  the  base  current  is  mainly  the 
diffused  hole  current  in  the  bulk  region  of  a  constant  doped  emitter.  Figure  2.10 
shows  the  hole  densities  in  the  emitter  region  with  different  interface  thick- 
nesses. It  can  be  seen  that  the  hole  densities  are  lower  in  the  entire  emitter 
region  for  a  thicker  interface  layer.  Though  there  are  holes  accumulated  in  the 
interface  side  of  the  emitter,  the  hole  densities  determining  the  base  current 
are  dominated  by  the  holes  injected  from  the  base-emitter  junction.  Thus,  the 
carrier  accumulation  and  depletion  at  the  edge  of  the  interface  is  not  an  impor- 
tant effect  for  the  reduction  of  base  and  collector  currents  at  high  biases. 

Another  effect  of  the  interface  insulator  layer  is  to  introduce  an 
additional  voltage  drop  across  the  insulator  due  to  the  charge  distribution  at 
both  sides  of  the  layer.  The  voltage  drop  can  be  equivalent  to  an  extra  emitter 
resistance.  Figure  2.11  shows  the  simulated  emitter  resistance  as  a  function  of 
Vbe  with  an  interface  thickness  of  10  A  .  Due  to  the  larger  carrier  tunneling 
probability  through  the  interface  at  higher  biases,  the  emitter  resistance 
decreases  as  Vbe  increases.  The  extra  voltage  drop  in  the  interface  causes  the 
reduction  of  the  voltage  across  the  base-emitter  junction.  The  total  voltage 
reduction  across  the  base-emitter  junction  AV  can  be  obtained  from  the  actual 
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Figure  2.11.  The  resistance  of  the  polyemitter  interface  from  theo  FLOODS 
simulations.  The  thickness  of  the  interface  is  10A  and  the 
emitter  area  is  100|im. 
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current  deviation  from  the  ideal  base  current  in  the  Gummel-plot  data  and  is 
given  by 

AV  =  RBIB  +  REIE  (2-18) 

At  low  current  injection,  all  the  curves  in  figure  2.9  are  close  together  because 
the  reduction  of  the  voltage  across  the  base-emitter  junction  is  not  significant. 
When  the  current  injection  level  becomes  high,  thicker  interface  with  higher 
emitter  resistance  leads  to  a  larger  AV  and  a  smaller  voltage  across  the  base- 
emitter  junction.  Thus,  the  hole  density  at  the  emitter  side  of  base-emitter 
junction  decreases  following  the  increase  of  the  interface  thickness  as  can  be 
seen  in  figure  2.10.  Since  AV  has  the  same  effect  to  base  and  collector  current 
of  a  specific  interface  thickness,  the  base  and  collector  currents  begin  to  deviate 
from  the  ideal  curves  at  the  same  bias. 

2.5  Device/Circuit  Tightly  Coupled  Simulator 
2.5.1  Motivation 

In  high  speed  circuit  design,  circuit  simulators  such  as  SPICE  are 
heavily  used.  One  problem  often  faced  is  the  inaccuracy  of  the  active  device 
models  in  the  simulators  which  either  require  time  intensive  parameter  extrac- 
tions or  nonphysical  parameter  fitting  steps  from  measurement  or  device 
simulations.  On  the  other  hand,  the  device  simulators  such  as  PISCES  provide 
more  accurate  predictions  of  a  single  active  device  with  predetermined  bound- 
ary conditions.  However,  the  application  of  the  device  simulations  to  the  circuit 
design  is  limited  since  it  is  difficult  to  predict  the  interaction  between  devices 
while  the  boundary  conditions  are  not  determined.  The  combination  of  device 
and  circuit  simulators  is  thus  preferred  in  this  respect. 
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For  modern  scaled  devices,  especially  the  high  speed  and  high  frequency 
devices,  the  resistances  and  capacitances  in  the  active  region  are  greatly 
reduced.  However,  the  parasitics  outside  the  active  region  are  not  scaled  down 
accordingly  and  become  important  in  high  frequency  operation.  To  properly 
predict  the  device  performances,  the  inclusion  of  all  the  structures  is  necessary. 
Another  advantage  of  the  combination  of  device  and  circuit  simulators  is  to  cut 
down  the  number  of  grid  points  in  the  simulated  device  structure.  Very  often  a 
structure  contains  purely  resistive  or  nondistributed  capacitive  sections  which 
can  be  represented  by  lumped  resistances  or  capacitances.  The  replacement  of 
these  massive  resistive  or  capacitive  structures  with  lumped  circuit  compo- 
nents greatly  reduces  the  number  of  grid  points  which  in  term  saves  cpu  time 
and  memory. 

The  coupled  device/circuit  simulators  have  been  applied  to  the  stead- 
state  and  transient  simulations  of  BJT  and  MOS  circuits  [R0I88,  Oka88]. 
However,  the  capability  of  a.c.  simulations,  one  of  the  most  important  applica- 
tions, is  not  available  so  far  in  this  type  of  simulator.  In  this  subsection,  the 
implementation  of  the  coupled  device/circuit  simulator  and  its  applications  to 
d.c.  and  a.c.  simulations  are  presented. 

2.5.2  Algorithms 

There  are  two  general  methods  to  implement  the  device/circuit  coupled 
simulator.  "Loosely-coupled"  approach  [Oka88,  May92,  Kiz93]  calculates  the 
terminal  informations  (current  I,  voltage  V)  of  the  devices  from  the  device 
simulator  and  the  derivatives  (dV/dl)  numerically.  The  circuit  simulator  takes 
the  informations  from  device  simulations  instead  of  an  internal  compact  model 
to  compute  the  node  voltages  and  branch  currents  in  the  network.  Several 
iterations  between  device   and  circuit  simulators   are  required  to  get  the 


33 


V 


Nl 


Circuit  Simulator  (V^,  Ib) 


V 


N2 


V- 


N3 


V, 


el 


Lel 


-r~£ 


V, 


e2 


Le2 


V, 


e3 


le3 


Electrode 
i  Mesh 


Interface 

Mesh 


Semiconductor 
Mesh 


Device  Simulator 


Circuit  /  Device  Tightly  Coupled  Simulator 


Figure  2.12.  The  diagram  of  the  coupling  algorithm  between  circuit  and  device 
simulators.  The  solution  variables  in  device  are  \\f,  n,  p,.,  and  the 
solution  variables  in  circuit  are  V^  and  Ig.  The  coupling  is 
through  the  solution  variables  Vei  and  Iei,  where  i  =  1,  2,... 
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Figure  2.13.  The  Jacobian  matrix  of  the  coupled  device/circuit  simulator.  The 
sub-matrices  S  and  I  are  assembled  from  the  semiconductor  and 
interface  boundary  equations  in  the  device  simulator, 
respectively,  and  the  sub-matrix  C  is  from  the  modified  nodal 
equations  in  circuit  simulator.  An  additional  sub-matrix  E  is 
assembled  to  couple  the  solutions  between  circuit  and  device. 
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solutions.  This  method  is  easy  to  implement  but  iterations  between  circuit  and 
device  simulations  are  time  consuming  and  may  cause  divergence  if  the  initial 
guess  and  the  derivatives  of  the  device  terminal  characteristics  are  not  well 
estimated. 

A  "tightly-coupled"  method  [R0I88]  solves  circuit  and  device  equations 
simultaneously  using  fully  coupled  and  self-consistent  algorithms.  This 
method  requires  no  iteration  between  device  and  circuit  simulations  and  exhib- 
its better  convergence  than  loosely-coupled  method.  Due  to  the  features  of 
faster  convergence  and  stability  the  tightly-coupled  algorithm  is  used  in  our 
device/circuit  simulator. 

The  algorithm  of  the  tightly-coupled  device/circuit  simulator  is  illus- 
trated in  figure  2.12.  The  circuit  simulator  is  implemented  with  modified  nodal 
equations  (MNE)  and  the  solution  variables  are  the  nodal  voltages  V^  and 
branch  currents  Ig  of  voltage  sources.  The  basic  structure  of  the  device  simula- 
tor is  intact  except  the  governing  equations  in  the  interface  and  electrode 
meshes.  The  boundary  condition  of  Poisson  equation  (RE.)  in  the  interface  is 
given  by 

EF  +  Vei  =  0  (2-19) 

where  Vei  is  the  contact  voltage  which  is  a  solution  variable  in  the  electrode 
mesh  and  coupled  to  the  circuit  nodal  voltage  Vpji  by  the  equation  in  the 
electrode  mesh 

The  contact  current  Iei  is  another  solution  variable  in  the  electrode  mesh  which 
is  also  a  branch  current  of  circuit  node  Ni.  The  coupling  of  Ie(  with  device  is 
given  by  the  governing  equation  in  the  electrode  mesh 
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hi    ~   Li'C£JJni  +  Jpi+Jdis,i) 


(2-21) 


where  L^  is  the  contact  length  which  converts  the  current  density  of  2D  device 
simulator  (in  unit  of  A/|im)  to  current  that  flows  through  the  circuit  network  if 
the  uniform  current  densities  at  the  contacts  of  device  are  assumed.  The  result- 
ing Jacobian  matrix  of  the  coupled  simulator  and  the  governing  equations  used 
are  plotted  in  figure  2.13.  The  sub-matrix  S  represents  the  semiconductor  and 
I  is  the  interface  Jacobian  from  the  device  simulator  and  sub-matrix  C  is  the 
Jacobian  of  circuit  simulator.  The  link  of  these  two  systems  is  through  the 
variables  Vej  and  Iej  in  the  electrode  mesh  of  the  device  simulator.  Vej  and  Iej 
serve  as  the  terminal  solutions  and  boundary  conditions  of  the  device  as  well 
as  the  initial  values  and  solutions  of  the  circuit  nodes  simultaneously.  The 
strong  coupling  between  these  equations  and  solution  variables  ensures  the 
fast  convergence  of  the  algorithm. 

In  FLOODS,  the  a.c.  simulation  scheme  is  implemented  with  the  sinuso- 
idal steady-state  analysis  techniques    [Lau85].   The  nodal   equations   of  a 
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and  inductor  L  is  vh  -  vt  -  L  ■  -=—  ,  where  v^ 
and  vt  are  the  small-signal  nodal  voltages  at  each  side  of  C  or  L,  and  ic  and  i^ 
represent  the  small-signal  branch  currents  of  C  and  L,  respectively.  The  Stiff- 
ness matrix  expressions  for  the  sinusoidal  steady-state  solutions  of  C  andL  are 
given  respectively  by 
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After  assembling  the  global  a.c.  matrix,  the  a.c.  system  is 


[J  +  jd\ 


X  =  B 


(2-22) 


where  J  is  the  d.c.  Jacobian  matrix,  D  contains  Stiffness  matrices  of  time  deriv- 
ative terms  such  as  those  of  C  and  L,  X  is  the  a.c.  solution  vector,  and  B  is  the 
vector  contains  the  a.c.  voltage  and  current  sources  of  the  circuit.  By  splitting 
the  system  into  real  and  imaginary  part,  eq.  (2-22)  can  be  rewritten  using  only 
real  arithmetic  as 

(2-23) 

where  subscript  R  indicates  the  real  part  and  /  is  the  imaginary  part  of  the 
vector.  The  solutions  of  the  small-signal  analysis  can  be  obtained  by  solving  the 
matrix. 
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2.5.3  D.C.  Simulations 

A  long-stripe,  single-metal,  double-poly-base-contacts  BJT  (see  figure 
2.14)  is  used  for  the  studies.  This  device  is  designed  by  Harris  Semiconductor 
Inc.  for  radio  frequency  applications.  Usually,  only  the  active  region  of  the 
device,  as  shown  in  figure  2.15,  is  simulated.  Conventional  simulations  assume 
the  device  structure  is  symmetric  to  the  center  of  the  emitter  and  the  biases  at 
both  base  poly  contacts  are  the  same.  Thus,  only  half  of  the  BJT  structure  is 
required  in  the  simulations  which  can  get  the  same  solutions  as  the  whole 
active  region  simulations  but  the  cpu  time  is  about  60%  less.  The  validity  of  the 
assumption  will  be  discussed  in  later  sections. 

Figure  2.16  shows  the  simulation  diagram  of  the  device  with  lumped 
parasitics,  where  Rpc  represents  the  resistance  between  metal  and  poly  base 
contacts  and  is  determined  from  the  measured  poly  sheet  resistance,  Rc  is  the 
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Figure  2.14.  Top  view  of  a  long  stripe  double  poly  base  BJT  used  in  this  work.  The  arrow 
signs  indicate  the  low  and  high  resistive  current  paths  for  NPOLY  and 
FPOLY,  respectively. 
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Figure  2.15.  The  cross  section  diagram  of  a  single  metal,  double  poly  base 
contacts  BJT.  The  area  within  the  dashed  line  is  the  simulation 
region. 
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Figure  2.16.  Schematic  diagram  of  a  BJT  with  parasitics  used  in  the  symmetric 
simulations.  The  NPOLY  and  FPOLY  contacts  are  shorted.  The  total 
current  flowing  through  each  contact  of  the  active  device  IP  equals  the 
current  density  JP  times  the  actual  contact  length  Lactual.  i.e.  IP  =  Lactual  * 
JpThe  parasitic  values  are  for  a  device  with  50u.m  emitter  stripe. 
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collector  resistance  in  the  buried  layer  as  well  as  contact  and  estimated  from 
additional  simulations,  Ccb  takes  into  account  the  capacitance  between  base 
poly  and  collector  layer,  C^  is  the  capacitance  between  base  and  emitter  poly 
layers,  and  Ccs  and  Rsub  are  the  series  capacitance  and  resistance  connecting 
to  substrate.  These  capacitive  parasitics  are  obtained  from  separate  measure- 
ments and/or  simulations.  The  connected  resistance  Rpe  between  NPOLY  and 
FPOLY  is  neglected  which  results  in  the  symmetric  biases  at  the  poly  base 
contacts  as  assumed  in  the  conventional  simulations.  The  current  density  in 
each  contact  of  the  device  is  multiplied  by  the  contact  length  (=  2W)  to  repre- 
sent the  current  flowing  in  the  circuit  network. 

The  comparison  of  Gummel-plot  between  different  simulations  and 
measured  data  is  shown  in  figure  2.17.  In  intrinsic  simulation,  neither  parasit- 
ics nor  polyemitter  interface  model  is  included.  The  base  current  of  the 
intrinsic  simulation  is  basically  a  straight  line  until  the  current  crowding  effect 
in  the  intrinsic  base  region  occurs.  Though  it  properly  predicts  the  current  at 
low  biases,  intrinsic  simulation  overestimates  both  the  base  and  collector 
currents  at  high  level  current  injection.  When  taking  the  parasitics  into 
account,  the  base  and  collector  currents  deviate  from  the  intrinsic  curves  due 
to  the  additional  voltage  drop  in  Rpc  which  reduces  the  voltage  across  the  base- 
emitter  junction.  The  inclusion  of  parasitics  improves  the  prediction  but  a  large 
discrepancy  is  still  observed  in  the  simulation.  If  the  polyemitter  model  with  a 

o 

10  A  interface  is  also  included,  both  base  and  collector  currents  further  bend 
down  due  to  the  extra  voltage  drop  in  the  interface  and  are  in  good  agreement 
with  the  measured  data.  Thus,  both  the  external  base  resistance  and  polyemit- 
ter interface  are  involved  in  the  high  current  effect  and  need  to  be  included  in 
the  simulations  for  accurately  predicting  the  d.c.  characteristics  of  devices. 
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Figure  2.18.  The  s-parameters  comparisons  between  Measured  data,  intrinsic 
device  simulation,  and  simulation  with  parasitics  (Symmetric 
structure),  (a)  input  reflection  coefficient  Sn. 
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Figure  2.18.  (b)  Reverse  transmission  coefficient  S12. 
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Figure  2.18.  (c)  Forward  transmission  coefficient  S2i- 
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Figure  2.18.  (d)  Output  reflection  coefficient  S22- 
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2.5.4  A.C.  Simulations 

To  characterize  the  high  frequency  performance,  scattering  parameters 
(s-parameters)  are  used  because  they  are  easier  to  measure  at  high  frequencies 
than  other  kinds  of  parameters  [Sze81].  S-parameters  are  measured  using  a 
matched  termination.  For  example,  to  measure  the  S1;L  we  measure  the  input 
reflection  coefficient  when  the  output  port  is  terminated  with  an  impedance 
equal  to  the  characteristic  impedance  of  the  transmission  line.  Using  the 
matched  termination,  a  traveling  wave  incident  on  the  load  will  be  totally 
absorbed  and  no  energy  will  be  returned  to  the  output  port.  Using  the  s-param- 
eters to  measure  a  transistor  has  an  advantage  that  the  transistor  does  not 
oscillate.  In  contrast,  if  we  were  to  use  a  short-  or  open-circuit  test,  the  transis- 
tor could  become  unstable  [Gon84] .  However,  it  is  difficult  to  properly  describe 
the  transmission  line  with  finite  RLC  circuit  elements  in  the  circuit  simula- 
tions. Thus,  the  y-parameters  are  solved  in  our  device/circuit  simulations  and 
then  converted  equivalently  to  s-parameters  [Gon84] . 

The  structure  shown  in  figure  2.16  is  used  for  the  a.c.  simulations.  The 
capacitive  parasitics  C^  and  Ccs  are  adjusted  to  get  the  best  fit  of  the  measured 
data  due  to  the  difficult  capacitance  extraction.  The  resistive  parasitics  are 
intact  in  the  simulations  since  the  extraction  is  straight-forward.  The 
simulated  s-parameters  of  a  common-emitter  configuration  B  JT  are  plotted  in 
figure  2.18  together  with  the  measured  data.  It  is  obvious  that  none  of  the  s- 
parameters  is  well  predicted  by  the  intrinsic  simulations.  The  simulated 
maximum  available  power  gain  frequency  fmag  is  19.5GHz  which  is  much 
higher  than  10.3GHz  of  the  measurement.  With  the  coupled  device/circuit 
simulator,  the  parasitic  effects  can  be  taken  into  account  and  the  predicted  fmag 
is  11.1GHz.  The  accuracy  of  the  simulated  s-parameters  is  significantly 
improved  as  can  be  seen  in  figure  2.18. 
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The  sensitivity  of  the  simulations  to  the  parasitics  is  also  analyzed.  It  is 
found  that,  among  the  parasitics,  the  base  external  resistance  Rpc  (NPOLY 
resistance)  is  the  most  important  parameter  which  strongly  affects  all  s- 
parameters.  Since  the  input  reflection  coefficient  Sjj  depends  only  on  the  input 
impedance  of  the  device  which  is  the  series  of  the  intrinsic  base  impedance  and 
Rpc,  thus,  Rpc  is  the  only  parasitic  controlling  S^.  Cc^  connects  the  base  and 
collector  electrodes  which  is  important  to  the  reverse  transmission  coefficient 
S}2  at  high  frequencies  due  to  the  direct  conduction  between  input  and  output 
ports.  Though  the  effective  capacitance  of  Cc^  in  the  input  port  is  amplified  due 
to  the  Miller  effect,  Cc^  has  little  effect  on  the  forward  transmission  coefficient 
S2i  and  Sj2  since  the  Miller  capacitance  of  Ccj,  is  parallel  with  Cn,  the  sum  of 
the  diffusion  capacitance  in  base-emitter  junction  and  the  Miller  capacitance 
in  base-collector  junction,  which  predominates  the  input  capacitance.  Except 
Cpe,  the  output  reflection  coefficient  S22  is  sensitive  to  all  the  parasitics, 
especially  the  substrate  parasitics  Ccs  and  Rsub  as  well  as  the  base  resistance 
Rpc.  Due  to  the  parallel  with  the  much  larger  capacitance  CK,  Cpe  shows  no 
effect  on  any  s-parameter. 

The  number  of  node  in  the  simulated  device  mesh  is  4348  which  is  much 
larger  than  the  number  required  by  the  circuit  simulator.  Therefore,  the 
increase  of  the  running  time  in  the  coupled  device/circuit  simulations  is  indis- 
tinguishable compared  with  the  intrinsic  device  simulations.  Since  the  circuit 
elements  are  all  linear  in  the  simulations,  the  convergence  rates  also  are  the 
same  in  both  simulations. 

2.6  Summary 
Several  new  models  and  functions  have  been  implemented  in  FLOODS: 
The  single  node  self-heating  model  has  successfully  reduced  the  full  thermody- 


49 

namic  simulation  time  while  a  good  accuracy  was  achieved.  The  good 
convergence  ability  of  the  current  boundary  condition  has  also  been  demon- 
strated. The  implementation  of  the  poly  emitter  interface  model  has  provided  a 
better  description  and  boundary  condition  for  the  polyemitter  BJT  and  the 
improvement  of  the  Gummel-plot  data  has  been  obtained.  A  circuit  simulator 
has  been  developed  and  tightly  coupled  to  FLOODS.  By  taking  into  account  the 
proper  interface  and  parasitic  effects,  good  agreements  have  been  achieved 
between  simulations  and  d.c.  as  well  as  a.c.  measurements. 


CHAPTER  3 

THREE  DIMENSIONAL  BASE  DISTRIBUTED  EFFECTS 

OF  STRIPE  BJTS--BASE  RESISTANCE  AT  DC 


3.1  Introduction 

Base  resistance  is  important  for  characterizing  the  large-  and  small- 
signal  behavior  [Tan79],  switching  characteristics  [Spi657  Gau79]  and  noise 
performance  [Zie58,  Auf96]  of  high  speed  bipolar  transistors.  The  accurate 
determination  of  the  base  resistance  is  essential  for  B  JT  design  and  modeling. 
It  is  known  that  the  base  resistance  decreases  as  the  operating  current 
increases  due  to  base  conductivity  modulation,  base  push-out  in  high-level 
injection,  and  current  crowding  effects  [Yua88,  Sat90,Fus92].  The  crowding  is 
the  tendency  for  current  to  flow  at  the  edge  of  the  intrinsic  base  region  at  high 
current  injection  levels  due  to  a  voltage  drop  along  the  base  current  path  in  the 
intrinsic  base  region.  The  current  crowding  reduces  the  effective  base  current 
path  and  decreases  the  base  resistance.  Accurate  modeling  of  base  resistance 
is  complicated  due  to  its  distributed  nature  and  operating  point  dependence. 

Several  papers  have  discussed  modeling  and  computation  of  the  base 
resistance  based  on  the  debiasing  and  current  crowding  effect  in  the  two- 
dimensional  (2D)  intrinsic  base  region  of  BJTs  [Yua88,  Rhe95,  Gho65,  Ver91, 
Cho90]  with  the  assumption  that  the  voltage  and  current  distribution  in  the 
third  dimension  are  uniform.  Schroter  [Sch91]  proposed  a  purely  geometry- 
dependent  model  to  compute  the  external  base  resistance.  However,  in  contem- 
porary high  speed  double  polysilicon  BJTs,  a  long  stripe  (see  Fig.  2.14)  or  long 
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interdigitated  base  structure  is  often  used  to  obtain  high  current  drive  and  low 
base  resistance.  In  this  situation  similar  debiasing  effect  occurs  in  the  polysil- 
icon  region  as  in  the  intrinsic  base  region  which  causes  nonuniform  voltage  and 
current  distribution  in  the  third  dimension  of  the  device.  The  voltage  and 
current  distribution  characteristics  vary  at  different  base-emitter  biases 
resulting  in  the  modulation  of  base  resistance.  To  date,  the  debiasing  effect  in 
the  base  polysilicon  region  and  the  three  dimensional  (3D)  nonuniform  current 
distribution  in  the  base  along  with  its  relationship  to  the  base  resistance  at 
different  biases  has  not  been  studied. 

The  purposes  of  this  chapter  are  to  model  the  d.c.  distributed  effect  in 
both  the  intrinsic  and  base  polysilicon  regions  and  to  develop  a  method  for 
computing  the  base  resistance  from  the  quasi-3D  model.  Using  the  results, 
effects  of  scaling  and  silicide  technologies  on  the  base  resistance  and  device 
design  consideration  are  also  examined. 

3.2  Asymmetric  Distributed  Model  in  the  Active  Base  Region 

3.2.1  Modeled  device  structure 

Figure  2.14  shows  a  typical  stripe  BJT  with  single  metal,  double  polysil- 
icon base  contact  structure  which  will  be  used  to  study  the  3D  voltage  and 
current  distributed  effect.  The  base  current  in  the  polysilicon  region  can  be 
divided  into  two  major  components.  One  flows  from  the  metal  base  contact  and 
injects  directly  into  the  intrinsic  base  region  at  the  NPOLY  (Near  POLY)  side. 
The  other  component  flows  through  the  FPOLY  (Far  POLY)  before  it  is  injected 
into  intrinsic  base  region.  Although  the  base  current  density  flowing  from 
metal  contact  to  the  edge  of  the  NPOLY  (lineN  in  Fig.  2.14)  might  be  different, 
the  voltage  along  lineN  can  be  assumed  to  be  constant  (=  VN )  due  to  the  small 
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resistance  and  negligible  ohmic  voltage  drop  difference  in  the  short  and  wide 
current  path.  The  base  current  flowing  through  the  FPOLY  experiences  larger 
resistance  due  to  a  longer  and  narrower  current  path  compared  with  current  in 
the  NPOLY,  which  results  in  a  larger  debiasing  effect  in  the  FPOLY.  The 
voltage  varies  along  the  FPOLY  and  is  expected  to  be  lower  than  that  of  the 
NPOLY. 

To  simplify  the  problem,  we  assume  the  base  current  injects  from  the 
base  polysilicon  contact  to  the  intrinsic  base  region  boundary  through  a  lumped 
extrinsic  base  resistance  R^ext  without  any  current  loss  (see  Fig.  3.1).  The  base 
current  in  the  intrinsic  base  region  is  also  assumed  to  flow  only  in  the  x-direc- 
tion  and  has  no  current  component  in  the  y-direction.  If  the  stripe  is  long,  the 
base  current  contributed  from  the  short  edges  of  the  emitter  is  negligible.  The 
base  current  in  the  FPOLY  can  then  be  assumed  to  flow  only  in  y-direction  and 
conduct  from  the  NPOLY  through  a  lumped  resistance  Rpe  at  the  edges  of 
polysilicon  base  handle  as  shown  in  Fig.  2.14.  Thus,  the  device  can  be  treated 
as  the  integration  of  numerous  slices  of  2D  structures  in  the  y-direction. 

3.2.2  Intrinsic  Model  Development 

Based  on  the  assumptions  and  derivation  of  Hauser  [Hau64] ,  the 
governing  equation  of  the  base  current  in  the  intrinsic  base  region  Jb(x)  may 
be  written  as  a  second-order  differential  equation: 


1  Jb(x)  +  ^Jb(x)4-Jb(x)  =  0  (3-1) 


dx2  V0 


KT 

where  pb  is  the  intrinsic  base  resistivity,  V 0  =  is  the  thermal  voltage,  and 

x  is  the  coordinate  in  the  intrinsic  base  region.  The  solution  of  this  equation  is 
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Jb(x)  =  A  •  tanJA  ■  ^- ■  b[\  - 1|  (3-2) 


where  A  and  £  are  arbitrary  constants  of  integration.  The  constants  can  be 
evaluated  from  the  boundary  conditions  of  the  device.  There  is  a  point  I  in  the 
intrinsic  base  at  which  the  potential  is  a  minimum  and  the  base  current  is  zero, 
i.e.  Jb(l)  =  0 .  This  gives  B  =  I.  Equation  (3-2)  may  then  written  as 

2VQ 
Jb(x)  =  — -  •  Z •  tan(Z •  (I- x))  (3-3) 

Pfo 

where 

Z  -  %-  (3-4) 

The  voltage  in  the  intrinsic  base  is  then  given  by 

X 

V(x)  =  V(Q)-jjb(x)-pbdx  =  y(0)-2Voln[^COS^s((^")X)))  (3-5) 

o 

In  case  of  a  symmetric  bias,  i.e.  V(0)  =  V(L),  where  L  is  the  intrinsic  base  length 
beneath  emitter,  eq.  (3-5)  gives  I  =  L/2.  Substituting  this  into  eq.  (3-3),  we 
obtain 

Jb(L)  =  — -Z-  tan  Z-|    =  -Jb(0)  (3-6) 

Then  Z  is  obtained  from  the  solution  of  the  equation. 
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NPOLY 


FPOLY 


-      Path  to  extract  potential  in  device  simulation 
Base-emitter  metallurgical  junction  in  the  device 


Figure  3.1.  Schematic  diagram  of  the  active  region  of  a  double  polysilicon  base 
BJT.  Shaded  regions  are  the  polysilicon  emitter  and  base.  A  lumped 
extrinsic  base  resistance  Rbext  connects  the  polysilicon  contact  and 
intrinsic  base  region.  All  the  base  current  injected  into  the  intrinsic 
base  area  flows  through  Rbext.  The  dashed  line  indicates  the  path 
used  to  extract  potential  from  FLOODS  simulations. 
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Combining  eq.  (3-3)  and  (3-5)  at  x=L,  the  relation  between  boundary  current 
and  voltage  is  given  by 


If  |V(0)-V(L)|  «2V0,then 


^i^uy^i^u  o-io) 


The  validity  of  this  assumption  will  be  verified  in  section  3.3.2.  With  this 
approximation,  eq.  (3-8)  and  (3-9)  become 

Jbm  •  it^im)V(0)~V(L))+2-if  tai# )  (3"U) 

jb(L)  -  \^mf(h)  - vm + 1J£-  tan(f ))       (3-12) 

and  the  total  current  is 

Jbt  =  — —  ■  tanf^O  (3"13) 

9b  V  2  J 

3.2.3  Model  with  Extrinsic  Base  Resistance 

Assuming  all  the  base  current  in  the  intrinsic  base  region  boundary 
flows  through  a  lumped  extrinsic  base  resistance  Rbext  from  the  polysilicon 
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base  contact  without  any  loss,  as  shown  in  Fig.  3.1,  the  I-V  relation  is  then 
given  by 

VN  =  V(0)+Jb(0)-Rbext  (3-14) 

VF  =  V(L)-Jb(L)-Rbext  (3-15) 

where  VN  and  VF  are  voltages  at  NPOLY  and  FPOLY,  respectively.  Substitut- 
ing eq.  (3-11)  to  (3-13)  into  (3-14)  and  (3-15),  gives  the  terminal  current  density 
equations  as  a  function  of  terminal  voltage  VF(y). 

JbN(y)  =  Jb(0)  =  C,  •  VF(y)  +  C2  (3-16) 

JbF(y)  =  -Jb(L)  =  -C,  •  VF(y)  +  C3  (3-17) 

where  y  is  the  coordinate  along  the  polysilicon  base  region, 

Z 


pb  •  sin(ZL)^ 
Z 


C{  = — — j- —  (3-18) 


**«<    Kpysin(ZL) 


c    _  Rb^Jbt  -  (P6  ■  si!(ZL))+  2^f  "  ^(f )+  Pb  ■  sm(ZL) '  V* 


(3-19) 


and 


c        Rb'"Jbt  -  U  -  ^(ZL))+  '-^f  •  taP(f )-  p,  -  sin(ZL)  -  V» 

l+2RbA-^lkzL) 


(3-20) 
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Though  the  base  current  density  flowing  from  metal  contact  to  the  edge  of 
NPOLY  (lineN  in  Fig.  2.14)  might  be  different,  the  voltage  along  lineN  is 
assumed  to  be  constant  (=  VN )  due  to  the  small  resistance  and  ohmic  voltage 
drop  in  the  short  and  wide  current  path.  So  if  we  assume  pb  and  Rbext  are 
independent  of  current  crowding  at  each  bias  and  the  total  base  current  Jbt 
remains  unchanged,  then  C{,  C2,  and  C3  become  constants.  Thus,  the  base 
current  densities  injected  from  NPOLY  JbN(y)  and  FPOLY  JbF(y)  into  the 
active  base  region  are  a  linear  function  of  VF(y) ,  which  is  determined  by  the 
debiasing  effect  along  the  FPOLY. 

.3.2.4  Verification  and  Discussion 

A  real  single  metal  base  contact  (SMB)  npn  BJT  with  an  emitter  stripe 
length  of  50|im  is  used  for  this  study.  The  layout  and  cross-section  of  the  intrin- 
sic device  structures  are  illustrated  in  Fig.  2.14  and  2.15,  respectively.  The  2D 
doping  profile  was  determined  using  the  process  simulator  Athena,  and  the 
polysilicon  sheet  resistance  was  obtained  from  independent  measurements. 
The  other  resistances  and  boundary  condition  J^(0)  required  for  the  model  are 
determined  from  2D  device  simulations  using  FLOODS  [15],  [16].  Due  to  the 
Early  effect  and  base  conductivity  modulation,  the  base  width  W5  and  intrinsic 
base  resistivity  pb  vary  depending  on  the  current  injection  level  and  need  to  be 
extracted  from  simulations  at  each  bias.  The  required  parameters  for  the 
model  of  the  transistor  are  summarized  in  Table  3.1. 

Figure  3.2  compares  the  normalized  voltage  distribution  along  the 
intrinsic  base  region  between  the  model  and  FLOODS  simulations.  In  all  cases, 
VN  is  fixed  at  0.8V  and  VFis  swept  from  0.8V  to  0.7994V  Both  the  model  and 
simulations  show  similar  debiasing  effects  across  the  intrinsic  base  region.  The 
discrepancy  between  simulation  results  and  the  model  is  due  to  the  fact  that 
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Table  3.1  The  device  parameters  required  for  model 


parameters 

value 

Emitter  width  (Intrinsic  base 
length),  L 

1pm 

Intrinsic  base  width,  Wb 
@Vbe=0.6V,  Vcb=2V 

0.22pm 

Unit  length  extrinsic  base 
resistance,  Rbext 

705.3fi.pm 

Intrinsic  base  resistivity,  pb 
@Vbe=0.6V,  Vcb=2V 

L2Ka.jj.rn 

Polysilicon  sheet  resistivity,  pp 

16.7Q.pm 

Polysilicon  thickness 

0.288pm 

FPOLY  width 

6pm 

FPOLY  cross  section  area,  Ap 

1.73pm2 

FPOLY  connecting  resistor,  Rpe 

12.57Q 

Resistance  from  lineN  to  metal, 
Rpc 

13.74Q 

Emitter  stripe  length,  2 W 

50pm 
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Figure  3.2.  The  DC  potential  profile  in  the  intrinsic  base  region  at  different  VFs 
(X=1.0  |im).  In  all  cases,  VNs  (X=0|im)  are  fixed  at  0.8V.  The  symbols 
indicate  the  device  simulation  data  and  lines  are  the  modeling 
results. 


60 

the  path  used  to  extract  the  potential  profiles  in  the  simulations,  as  indicated 
by  the  dashed  line  beneath  the  emitter  in  Fig.  3.1,  does  not  correspond  to  the 
actual  base  current  path  assumed  in  the  model.  Note,  at  the  minimum  voltage 
point  I  of  the  potential  profile,  where  J^(x)  is  zero,  the  base  current  injected 
from  base  to  emitter  on  the  left  hand  side  of/  is  contributed  by  the  left  terminal 
contact  (NPOLY),  and  the  current  injected  on  the  right  hand  side  of  /  comes 
from  the  right  terminal  contact  (FPOLY).  The  point  /  acts  as  a  current  divider 
in  the  intrinsic  base  region.  When  the  bias  is  symmetric,  I  is  located  at  the 
center  of  the  intrinsic  base  and  the  current  is  equally  divided  between  NPOLY 
and  FPOLY.  As  Vp  decreases,  I  moves  toward  FPOLY  and  forces  part  of  the 
current  originally  flowing  through  the  FPOLY  at  symmetric  bias  to  flow 
through  the  NPOLY.  Consequently,  the  current  J^  from  the  NPOLY  increases 
corresponding  to  the  decrease  of  current  J^  through  the  FPOLY. 

The  leverage  behaviors  of  terminal  currents  at  asymmetric  biases  can  be 
observed  from  Fig.  3.3.  The  FPOLY  current  JbF  decreases  and  the  NPOLY 
current  J^n  increases  if  Vp  is  smaller  than  V}y.  In  the  case  of  V}y=0.8V,  the  base 
current  injected  from  base  to  emitter  is  approximately  constant  everywhere  in 
the  intrinsic  base  region  since  the  voltage  variation  in  the  FPOLY  and  in  the 
intrinsic  base  region  are  small.  Thus,  the  reduction  of  J^p  is  compensated  by 
the  increase  of  J^,  and  the  total  current  density  is  approximately  constant  at 
different  Vps  as  shown  in  Fig.  3.3(a).  As  the  bias  increases,  a  more  significant 
ohmic  voltage  drop  takes  place  along  the  FPOLY  which  reduces  the  voltage 
across  the  base-emitter  junction  and  decreases  the  total  base  current  density. 
Figure  3.3(b)  shows  the  J-VF  relation  while  VN  is  biased  at  0.92V.  The  total 
current  reduces  by  6.7%  when  Vp  is  6mV  lower  than  VN.  Under  these  high 
injection  conditions,  the  coefficient  Cj  in  Eq.  (3-16)  is  slightly  different  from  Cj 
in  Eq.  (3-17).  For  better  accuracy  in  the  quasi-3D  model,  CjS  in  Eq.  (3-16)  and 
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Figure  3.3.  Simulated  and  modeled  polysilicon  base  terminal  current  as  a 
function  of  VF  at  (a)  VN=  0.8V,  and  (b)  VN=  0.92V.  In  (b),  C1  in  (3.16) 
and  (3.17)  is  adjusted  to  meet  the  decreasing  total  current. 
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(3-17)  are  corrected  according  to  the  simulated  slopes  of  J^  and  J^F  at  high 
biases  as  illustrated  in  Fig.  3.3.  Though  it  requires  a  correction  in  the  coeffi- 
cients, the  linear  relation  between  the  terminal  currents  and  VF,  i.e.  Eq.(3-16) 
and  (3-17),  is  valid  until  the  voltage  difference  approaches  2V0  which  is  at  a 
bias  exceeding  the  range  of  interest. 

3.3  Quasi  Three  Dimensional  DC  Distributed  Model 

3.3.1  Model  Development 

The  following  assumptions  have  been  made  while  deriving  the  expres- 
sions for  current  and  voltage  distribution  in  the  polysilicon  base  regions:  (i)  the 
resistivity  of  polysilicon  base  p  is  constant,  (ii)  the  current  density  in  the 
FPOLY  is  uniformly  distributed  in  the  direction  perpendicular  to  the  current 
flow  (y-direction),  and  (iii)  the  device  is  symmetric  in  y-direction. 

Let  Jpp(y)  be  the  current  density  of  transverse  base  current  in  the 
FPOLY  and  be  related  to  the  terminal  current  JbF(y)  by  the  equation 

w 
JpF(y)  =  JpF(W)  +  j  —j^-dy  (3-21) 

y 

where  Ap  is  the  cross  section  area  of  the  polysilicon  base,  and  W  is  half  of  the 
emitter  stripe  length.  Since  the  device  is  symmetric  to  the  y-direction,  we  have 
JpF(W)  =  0 .  The  voltage  in  the  polysilicon  base  VF(y)  is  given  by 

y 
VF(y)  =  VF(0)-jjpF(y)ppdy  (3-22) 

o 

Substituting  eq.  (3-17)  and  (3-22)  into  (3-21)  and  differentiating  twice,  the 
result  is 
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,2 

JpF{y)-CAJpF{y)  =  0  (3-23) 


dy 


where 


Ap 


C4-J— 7^1  (3-24) 


Applying  the  boundary  condition  JpF(W)  =  0  to  the  general  solution  of  eq.  (3- 
23),  we  have 

JpF(y)  =  c5  ■  sinh(C4  •  (y  -  W))  (3-25) 

where  C5  is  a  constant  and  can  be  expressed  as 

JpF(0) 
°5  =  "sinh(C4W)  (3"26) 

If  eq.  (3-25)  is  used  in  (3-22)  and  the  integration  performed,  one  obtains 

VF(y)  =  C6-Cr  cosh(C4  -(y-W))  (3-27) 

where 

C7  =  4^  (3-28) 

and 

C6  =  V>(0)  +  C7cosh(C4  •  W)  (3-29) 

If  the  3D  nonuniform  current  flow  at  the  short  edge  of  the  emitter  stripe  is 
neglected,  and  we  assume  all  the  current  in  the  FPOLY  flows  through  the  resis- 
tance Rpe  at  the  edge  of  base  polysilicon,  the  edge  voltage  VF(0)  of  the  FPOLY 
and  Vjy  are  related  by 

VF(0)  =  VN-(Ap-Rpe)JpF(0)  (3-30) 
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Substituting  eq.  (3-26),  (3-27)  and  (3-30)  into  the  combination  of  (3-17)  and  (3- 
21),  and  the  FPOLY  boundary  current  is  given  by 

J^(0)  "  ApiC^WRpe)  ■ (C' - C< V">  (3-3 " 

which  is  only  a  function  of  V/y  where 

C8  =  l  +  -ir£--irz_coth(C4W).  (3-32) 

ApC4     APL* 

Once  Jpp(0)  is  obtained,  the  voltage  Vp(y)  and  base  current  density  Jpp(y)  in 
the  FPOLY  can  be  solved,  which  in  turn  gives  the  voltage  and  current  distribu- 
tion in  the  whole  area  of  the  intrinsic  base  region. 


3.3.2  Results  and  Discussion 

The  modeled  ohmic  voltage  drop  along  the  FPOLY  at  three  different 
biases  is  shown  in  Fig.3.4.  At  Vbe=0.8V,  the  model  shows  a  small  debiasing 
effect  along  the  FPOLY  and  the  voltage  variation  is  less  than  0.3mV.  When 
Vbe=lV,  the  ohmic  voltage  drop  in  the  FPOLY  is  larger  than  18mV  due  to  the 
larger  current  flow  in  the  FPOLY.  However,  even  at  this  condition, 
AV"  «  2VTH  ~  52mV ,  the  linear  approximation  in  deriving  Eq.(3-ll)  and  (3-12) 
is  still  valid.  The  emitter  current  density  injected  through  the  base-emitter 
junction  may  change  dramatically  in  the  presence  of  the  debiasing  effect 
because  of  its  exponential  dependence  on  the  base-emitter  junction  voltage. 
Figure  3.5  shows  the  modeled  emitter  current  density  Je  beneath  the  emitter 
area  when  Vbe  =  1.0V.  Je  is  constant  along  the  intrinsic  base  edge  near  NPOLY 
because  of  the  constant  bias  at  the  NPOLY  whereas  Je  decreases  toward  the 
center  of  the  stripe  due  to  the  ohmic  voltage  drop  along  the  FPOLY  in  the  y- 
direction  and  in  the  intrinsic  base  region  along  the  x-direction.  Thus,  the 
emitter  current  crowds  at  the  edges  of  the  stripe.  The  total  emitter  current 


65 


0 


<i 


■10 


■15 


-20 


0 


5  10  15  20 

Y  axis  along  FPOLY 


25 


Figure  3.4.  Modeled  DC  voltage  variation  in  FPOLY  at  three  different  biases. 
The  emitter  stripe  length  is  50jim  and  the  width  of  FPOLY  is  6|im 
with  an  average  sheet  resistance  of  58QJsq.  The  voltage  is  assumed 
to  be  symmetric  around  the  center  of  FPOLY. 
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Figure  3.5.  Plot  of  the  two-dimensional  distribution  of  emitter  current 
calculated  by  the  model.  The  intrinsic  base  sheet  resistance  is 
4.1KO/sq. 
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computed  from  the  quasi-3D  model  is  smaller  than  that  of  the  2D  model 
[Hau64]  due  to  the  additional  debiasing  effect  in  the  y-direction.  The  base 
current  density  is  proportional  to  the  potential  gradient  in  the  intrinsic  base 
region  of  xy-plane  if  the  conductivity  is  constant. 

Although  the  debiasing  in  the  FPOLY  at  low  biases  is  very  small,  it  still 
cause  the  voltage  in  the  FPOLY  to  be  slightly  lower  than  that  in  the  NPOLY. 
The  modeled  terminal  base  current  densities  injected  from  the  polysilicon  to 
the  active  base  region  along  the  y-direction  are  shown  in  Fig.  3.6.  The  current, 
JbN,  injected  from  the  NPOLY  is  always  larger  than  the  current,  JbF>  from  the 
FPOLY  since  the  bias  at  the  NPOLY  is  higher  than  the  bias  at  the  FPOLY.  The 
difference  between  JbN  and  JbF  is  larger  at  the  center  than  at  the  edge  of  the 
stripe  due  to  the  debiasing  along  the  FPOLY.  Note  that  the  decreasing  current 
in  the  FPOLY  is  compensated  by  the  increasing  current  in  the  NPOLY  if  VN  is 
less  than  0.9V.  The  total  current  is  about  the  same  as  that  without  considering 
the  debiasing  effect  in  the  FPOLY.  However,  when  VN  is  larger  than  0.9V,  the 
current  decrease  in  the  FPOLY  is  larger  than  the  increase  in  the  NPOLY 
because  of  the  significant  ohmic  voltage  drop  in  the  FPOLY.  The  total  current 
is  less  than  that  without  taking  into  account  the  debiasing  effect  in  the  FPOLY. 

3.3.3  Model  Verification/Application 

For  a  symmetric  double  polysilicon,  double  metal  base  contact  structure, 
the  base  current  contributed  from  each  side  of  the  polysilicon  base  should  be 
the  same.  However,  for  the  single  metal  base  contact  device  as  illustrated  in 
Fig.  2.14,  the  current  contributed  by  the  FPOLY  is  less  than  that  from  the 
NPOLY  due  to  the  asymmetric  bias  as  has  been  discussed  above.  Figure  3.7 
shows  the  ratio  of  base  current  contributed  from  the  FPOLY  and  NPOLY  as  a 
function  of  bias.  Note  that  even  at  low  biases,  the  ratio  is  less  than  40%  and 
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Figure  3.6.  Modeled  DC  base  terminal  current  density  injected  into  active  base 
region  vs.  position  along  base  polysilicon  at  different  biases.  The 
stripe  length  is  50\im  and  voltage  in  the  NPOLY  is  assumed  to  be 
uniform.  The  current  densities  JbN  and  Jy?  are  normalized  to  that 
aty=0.0(j.m. 
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Figure  3.7.  The  ratio  of  base  current  contributed  by  FPOLY,  IPF,  to  the  total 
current,  Ipp+Ippj,  as  a  function  of  bias  calculated  by  the  quasi-3D 
model,  where  IPN  is  the  current  contributed  by  the  NPOLY. 
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falls  off  slightly  as  the  bias  increases.  For  biases  larger  than  0.9V,  the  ratio 
decreases  very  fast  due  to  the  significant  debiasing  effect  in  the  FPOLY. 

The  current  flowing  through  the  FPOLY  experiences  a  larger  base  resis- 
tance than  that  through  the  NPOLY.  Since  the  current  via  FPOLY  decreases 
and  the  current  associated  with  NPOLY  increases  at  higher  biases,  it  implies 
more  and  more  current  flows  through  the  low  resistive  path  as  the  bias 
increases.  Thus,  both  the  current  crowding  effect  in  the  intrinsic  base  region 
and  the  current  redistribution  in  the  polysilicon  base  region  contribute  to  the 
decrease  of  base  resistance  as  the  bias  increases. 

The  total  base  resistance  R^  can  be  calculated  using  the  dissipated 
power  method  [Dub94] 

*6  =  ^**  =  ZjD|«V  W*  '  <*»i  (3-33) 

where  i>;  is  the  volume  of  the  region  of  interest,  7j,  is  the  total  base  current,  Jj, 
is  the  current  density  flowing  through  v^,  and  VV&  is  the  gradient  of  potential 
in  the  direction  of  current  flow.  Figure  3.8  shows  the  division  of  regions  used  to 
compute  the  individual  power  dissipation  in  each  region.  The  total  power  dissi- 
pation, P^  is 

°6  =  Pbif  +  Pbin  +  Ppf  +  Pben  +  Pbef  +  Ppc  +  PpR  (3-34) 

where  Ppc  =  IbRPC  is  the  power  dissipated  in  NPOLY  through  a  lumped 
resistance  Rpc-  Ppe  =  1R  •  (Ap  JpF(0))  is  the  power  loss  in  the  polysilicon 
connection  between  NPOLY  and  FPOLY.  Pben  anc^  ^BEF  are  dissipated  power 
in  the  extrinsic  base  regions  due  to  lumped  resistances,  R{,ext,  connected  to 
NPOLY  and  FPOLY,  respectively.  Ppp  is  the  power  consumption  in  the  FPOLY, 
PBjp- and  PBIN  are  the  power  loss  in  the  intrinsic  base  regions  due  to  the  current 
flow  associated  with  FPOLY  and  NPOLY,  respectively.  Since  the  base  current 


71 


oxide 


Rbext  ^bin       ^bif     !   ^bext 


_X 


Figure  3.8.  The  diagram  of  the  long  stripe  BJT  structure  and  the  definition  of 
the  base  resistance  in  each  subregion.  The  lightly  shaded  area  is 
the  base  polysilicon  region.  Rpc,  Rpe,  and  R^ext  are  the  constant 
lumped  resistances.  Rnn,  Rfcf  and  RpF  are  functions  of  the  base 
current. 
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and  voltage  in  each  region  have  been  obtained  from  the  model,  the  total  dissi- 
pated power  and  base  resistance  can  be  solved  easily. 

The  measured  and  modeled  large-signal  base  resistances  are  plotted  in 
Fig.  3.9  as  a  function  of  bias.  The  Gummel-plot  method  [Fus95]  is  used  to 
measure  the  DC  base  resistance.  The  2D  model  curve  indicates  the  resistance 
calculated  with  NPOLY  and  FPOLY  shorted  together  which  results  in  the 
symmetric  voltage  and  current  in  both  sides  of  the  emitter  stripe.  Under  this 
condition,  Pben=pBEF,  PBIN=PBIF>  and  PPF  =  °-  Tne  base  resistance  of  the 
quasi-3D  model  is  seen  to  be  much  larger  and  more  accurate  than  that  of  the 
2D  model  in  the  bias  range  of  interest.  The  analysis  of  the  power  dissipation  in 
each  region  shows  the  difference  comes  from  two  facts,  (i)the  current  flowing 
through  NPOLY  in  a  quasi-3D  model  has  longer  current  paths  and  higher  effec- 
tive resistance  than  that  of  a  2D  model  due  to  the  movement  of  the  minimum 
voltage  point  toward  FPOLY.  (ii)  the  highly  resistive  current  path  of  the  FPOLY 
in  the  quasi-3D  model. 

As  the  bias  increases,  the  base  resistance  of  the  2D  model  mainly 
decreases  due  to  current  crowding  and  conductivity  modulation  in  the  intrinsic 
base  region.  In  addition  to  the  these  mechanisms,  the  current  redistribution 
from  high  to  low  resistive  paths  in  the  polysilicon  region  further  reduces  the 
base  resistance  in  the  quasi-3D  model  as  the  bias  increases.  The  variation  of 
base  resistance  computed  by  the  quasi-3D  model  is  thus  larger  and  more 
accurate  than  that  of  the  2D  model  as  can  be  seen  in  Fig.  3.9. 

The  2D  model  represents  a  symmetric  double  metal  base  contact  device 
structure  (DMB)  and  quasi-3D  model  describes  a  single  metal  base  contact 
structure  (SMB).  Figure  3.10  shows  the  modeled  base  resistance  plotted  as  a 
function  of  emitter  stripe  ratio  5R  with  three  different  FPOLY  widths,  X,  while 
Vbe  is  biased  at  0.9V.  As  expected,  the  base  resistance  of  DMB  decreases 
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inversely  proportional  to  the  increase  of  stripe  length.  The  Rb  of  SMB  always 
has  higher  values  than  that  of  DMB  and  the  device  with  a  narrower  X  results 
in  a  larger  Rb.  The  difference  of  Rb  between  these  two  structures  is  also 
included  in  Fig.  3.10.  It  is  shown  that  the  difference  is  small  for  the  shorter 
stripe  device  and  becomes  larger  as  the  stripe  length  increases  because  the 
distributed  effects  are  more  significant  at  a  longer  FPOLY.  Scaling  of  the 
FPOLY  width  further  enlarges  the  difference.  For  example,  if  9t  is  10,  Rb  of 
SMB  with  fc=l\im  is  26%  higher  whereas  Rb  of  SMB  with  h=2\im  is  23%  higher 
than  Rb  of  DMB. 

Silicided  polysilicon  technologies  are  often  utilized  in  contemporary  high 
speed  double  polysilicon  BJTs  to  reduce  the  base  resistance.  However,  this 
process  increases  the  complexity  of  process  and  cost.  Figure  3.11  shows  the 
difference  of  Rb  between  DMB  and  SMB  structures  computed  with  different 
polysilicon  sheet  resistances  and  SR  .  For  a  nonsilicided  polysilicon  with  a  sheet 
resistance  of  100  Q/D  ,  Rb  of  SMB  with  <R  =6  is  21%  higher  than  Rb  of  DMB. 
The  difference  increases  to  36%  if  SK=14.  If  the  sheet  resistance  is  reduced  to 
5 Q/D  ,  the  difference  decreases  but  is  still  significant  for  large  9t .  Therefore, 
the  design  of  SMB  structure  using  symmetric  device  simulations  without 
taking  into  account  the  distributed  effects  can  significantly  underestimate  the 
base  resistance.  The  results  also  suggest  the  Rb  of  DMB  device  is  always  much 
smaller  than  SMB  structure  unless  a  low  sheet  resistance  polysilicon  is  used. 
The  silicided  technology  is  thus  very  important  for  the  performance  of  SMB 
structure. 

3.4  Summary 
A  quasi-3D  model  describing  the  distributed  base  voltage  and  current  in 
both  the  intrinsic  and  polysilicon  base  regions  has  been  derived.  The  results 
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illustrate  the  nonuniform  characteristics  of  voltage  and  current  in  a  long  stripe 
transistor  that  cannot  be  properly  described  by  a  2D  model.  A  current  redistri- 
bution effect  in  the  polysilicon  base  regions  due  to  the  asymmetric  bias 
modulation  has  been  presented.  In  addition  to  the  current  crowding  effect  in 
the  intrinsic  base,  this  current  redistribution  effect  also  plays  a  key  role  in 
determining  the  bias  dependent  base  resistance.  The  base  resistance  has  been 
computed  from  the  solution  of  the  model  and  the  accuracy  has  been  verified  by 
the  measured  data.  The  same  concept  and  methodology  can  be  applied  to 
analyze  the  distributed  effects  in  the  long  interdigitated  structures. 
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Figure  3.9.  A  comparison  of  DC  base  resistances  obtained  from  measurements, 
2D,  and  quasi-3D  base  distributed  models.  The  emitter  stripe 
length  is  50)im  and  the  width  of  the  FPOLY  in  the  device  is  6(im. 
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Figure  3.10.  The  plot  of  the  modeled  base  resistance  and  difference  between 
SMB  and  DMB  structures  as  a  function  of  emitter  stripe  ratio  with 
different  FPOLY  widths.  The  sheet  resistance  of  FPOLY  is  58Q/sq 
and  FPOLY  width  is  indicated  in  the  parentheses. 
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Figure  3.11.  The  difference  of  base  resistance  between  SMB  and  DMB 
structures  as  a  function  of  emitter  stripe  ratio  with  different 
FPOLY  sheet  resistances.  The  width  of  FPOLY  is  l|im  and 
Vbe=0.9V. 


CHAPTER  4 

THREE  DIMENSIONAL  BASE  DISTRIBUTED  EFFECTS 

OF  STRIPE  BJTS-AC  EFFECTS  ON  INPUT  CHARACTERISTICS 


4.1  Introduction 

Design  considerations  for  a  microwave  amplifier  include  stability,  power 
gain,  bandwidth,  noise  and  d.c.  bias  requirements  [Gon84].  A  construction  of 
an  integrated  circuit  microwave  amplifier  usually  starts  with  a  proper  design 
of  the  transistor;  then  the  matching  networks  of  the  input  and  output  ports  are 
determined  for  a  particular  stability  and  gain  criteria.  An  unconditionally 
stable  transistor  will  not  oscillate  with  any  passive  termination.  On  the  other 
hand,  a  design  using  potentially  unstable  transistors  requires  careful  analysis 
so  that  the  matching  network  produces  a  stable  amplifier.  Thus,  the  accurate 
determination  of  the  transistor  characteristics  is  the  first  step  to  optimize  the 
performance  of  a  microwave  amplifier. 

The  input  impedance  of  a  transistor  in  the  common-emitter  configura- 
tion is  usually  difficult  to  predict  either  by  compact  models  or  2D  device 
simulations  due  to  the  3D  distributed  effect  in  the  base  region.  This  effect 
decreases  the  base  impedance  in  magnitude  with  increasing  frequency  as  a 
result  of  a  high-frequency  internal  biasing  effect.  Several  authors  [Gho65, 
Pri61,  Ver91]  have  modeled  the  base  impedance  with  a  uniform  transmission 
line  in  the  intrinsic  base  region  to  account  for  the  a.c.  distributed  effect  at  high 
frequency  operations.  However,  in  double  polysilicon  bipolar  transistors,  the 
distributed  characteristic  in  the  polysilicon  base  region  can  also  be  significant 
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at  high  frequencies  and  no  studies  of  this  effect  have  been  reported.  Neglecting 
this  effect  may  result  in  an  inaccurate  estimation  of  input  impedance  in  the 
transistor  and  a  poor  design  of  an  amplifier. 

In  this  chapter,  we  concentrate  on  the  distributed  effect  in  the  polysilicon 
base  region,  its  effect  on  the  input  characteristics,  and  the  accurate  and 
efficient  simulations  of  a  stripe  BJT  at  high  frequencies.  Since  the  distributed 
effect  is  a  3D  problem,  a  quasi-3D  simulation  scheme  is  proposed  to  character- 
ize the  input  behavior  of  the  device.  In  the  scheme,  FLOODS  is  combined  with 
an  analytical  model  to  deal  with  the  voltage  and  current  nonuniformity  in  the 
third  dimension  of  the  device  due  to  the  distributed  effects  in  the  base  polysil- 
icon region.  The  comparisons  between  quasi-3D,  conventional  2D  simulations 
and  the  measured  data  are  performed  and  the  results  are  discussed. 

4.2  A.C.  Distributed  Model  in  the  Active  Base  Region 
4.2.1  Model  Development 

Figure  4.1(a)  shows  a  distributed  network  representation  of  a  rectangu- 
lar active  base  region  if  the  emitter  and  collector  contacts  are  a.c.  grounded 
[Gho65].  The  diode  represents  the  junction  between  base  and  emitter  and  the 
capacitance  C^  accounts  for  the  base-collector  depletion  capacitance  under  the 
emitter  area  and  the  geometric  parasitic  capacitance  of  the  device.  At  a  specific 
d.c.  bias,  the  diode  can  be  linearized  by  a  parallel  differential  resistance  Rd  and 
a  capacitance  [Shr90]  if  the  series  inductance  is  negligible.  The  small-signal 
admittance  of  a  section  can  be  represented  by  Yd  =  Rd  +  jwCd  at  a  specific 
angular  frequency  w,  where  Cd  includes  the  diffusion  and  depletion  capacitance 
of  base-emitter  junction  and  C^.  The  complete  small-signal  network  is  plotted 
in  figure  4.1(b). 
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Figure  4.1.  Models  of  the  distributed  network  in  the  active  base  region,  (a) 
Large-signal  network,  (b)  Small-signal  network. 
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At  high  current  densities,  high-current  effects  like  Kirk  effect  [Kir62] 
cause  an  increase  of  the  transit  time  and  hence  a  reduction  of  the  cutoff 
frequency  fp  Because  of  this,  we  thus  limit  our  interests  in  the  current  range 
where  the  maximum  fy  is  achieved.  For  the  device  used  in  this  analysis,  the 
maximum  /y  is  achieved  at  Vbe  =  0.83V.  Under  this  condition,  d.c.  debiasing 
effects  in  the  intrinsic  base  as  well  as  polysilicon  base  regions  are  not  signifi- 
cant (see  figure  3.4  and  3.6)  and  the  current  injection  from  base  to  emitter  is 
approximately  constant  in  the  intrinsic  base  area.  Thus,  Y$  can  be  assumed  to 
be  uniform  in  the  intrinsic  base  region,  and  the  network  becomes  a  uniform  RC 
transmission  line. 

The  equations  governing  the  small-signal  voltage  vbi(x)  and  transverse 
current  density  jbi(x)  in  the  intrinsic  base  region  are  given  by 

j£Jbi{x)  =  -pw-  jbi(x)  (4-2) 

^JblM  =  -Yd-vbl(x)  («) 

where  p^j  is  the  average  sheet  resistance  of  the  base  layer.  Substituting  eq.  (4- 
3)  to  (4-2),  the  resulting  equation  is 

1  Jbtix)  =  YdPbljbl(x)  (4-4) 


dx 


The  solution  of  this  equation  is 

jbl(x)  =  C,  sinh(/yydPfotx)  +  C2cosh(jYdpblx)  (4-5) 

Then, 

vb^)  =  -  j^(Ci  cosh(  7*^*)  +  C2sinh(7l>^x))  (4-6) 
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where  Cj  and  C^  are  two  arbitrary  constants  to  be  determined  from  the  bound- 
ary conditions.  Solving  j^Cx)  and  vbi(x)  at  x=0,  it  gives 

Cl    =  -a/?UW(0)  (4"?) 

HPbi 

and 

C2  =  ,/M(0)  (4-8) 

Substituting  eq.  (4-7)  and  (4-8)  into  (4-5)  and  (4-6)  and  solving  at  x=L  gives  the 
small-signal  j- v  relations 

Jbitfi)  =  Cyvbi(0)-C4-vbi(L)  (4-9) 

and 

jbi(L)  =  C4-vbi(0)-C2-vbi(L)  (4-10) 

where 

C3  =  Jpcoth(VF^L)  (4-11) 

HPbi 

and 

IF 


d 


C4  =    Mc8ch(7FdP«L)  (4-12) 

If  the  current  densities  at  the  edges  of  intrinsic  base  region  flow  through  the 
lumped  extrinsic  base  resistance  Rbext  without  loss  and  equal  to  the  current 
densities  of  poly  contacts,  then  similar  equations  as  (3-14)  and  (3-15)  can  be 
used  to  transform  the  terminal  voltages  at  intrinsic  base  edges  to  that  at  base 
poly  contacts.  The  results  are: 
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jN(y)  =  JM(0)  =  C5-vN-C6-vF(y)  (4-13) 

JF(y)  =  -jbi(L)  =  C5  ■  vF(y)  -C6-vN  (4-14) 

where 

C3  +  C4C7Rbext 
C5  =  l+Rbext(C3  +  C4)  (4"15) 

C  a—C  AC  nRb  „„t 


and 


°i  =  l+Rbe3xt.(C3-C4)  (4"17) 


4.2.2  Verification  and  Discussion 

The  device  structure  shown  in  figure  2.14  and  2.15  is  used  for  the  a.c. 
analysis.  The  small-signal  voltage  vN  along  lineN  is  assumed  to  be  constant 
and  the  voltage  along  FPOLY  vp(y)  varies  due  to  the  debiasing  effect.  The 
small-signal  admittance  Yg  used  in  the  model  is  obtained  from  the  FLOODS 
simulations  of  two  p-n  junctions.  One  has  the  same  doping  profile  as  the 
modeled  BJT  in  the  base  and  emitter  regions  and  is  forward  biased  to  extract 
the  base-emitter  junction  admittance.  The  base-collector  depletion  capacitance 
is  extracted  from  the  reversed  biased  junction  whose  doping  profile  is  the  same 
as  the  modeled  BJT  in  the  base  and  collector  region.  The  geometric  parasitic 
capacitance  is  neglected  in  the  evaluation  of  the  admittance  Y^. 

The  normalized  amplitude  of  the  small-signal  base  current  densities 
injecting  from  the  intrinsic  base  region  to  emitter  are  plotted  in  figure  4.2  at 
different  frequencies  while  the  biases  are  symmetric.  The  current  crowding 
effect  is  significant  only  if  the  oscillation  frequency  is  higher  than  1GHz.  Figure 
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4.3  illustrates  the  same  plot  as  figure  4.2  but  with  asymmetric  biases.  It  can  be 
clearly  seen  that  the  position  of  minimum  voltage  /  moves  to  the  lower  bias 
terminal.  The  lower  the  frequency,  the  farther  the  point  I  moves.  For  the  cases 
the  operating  frequencies  are  1GHz  and  100MHz,  the  point  I  is  virtually  out  of 
the  intrinsic  base  region.  At  this  condition,  all  of  the  base  current  injecting 
from  base  to  emitter  is  supplied  by  the  higher  bias  terminal.  Part  of  the  current 
from  the  higher  bias  terminal  flows  to  the  lower  bias  terminal  via  the  intrinsic 
base  region.  The  current  direction  of  the  lower  bias  terminal  is  thus  changed. 
Figure  4.4  shows  the  small-signal  base  current  densities  at  the  poly 
terminals  as  a  function  of  vp(y)  oscillating  at  a  frequency  of  100MHz  when  Vjy 
is  biased  at  3mV.  If  vp(y)  is  slightly  smaller  than  v^,  the  real  part  of  the  low 
bias  terminal  current  density  Jp  decreases  and  the  reduced  amount  of  Jp  is 
compensated  by  the  increase  of  current  density  at  the  high  bias  terminal  J^ 
due  to  the  negligible  decrease  of  total  injection  current.  If  vp(y)  is  much  smaller 
than  vjy,  the  current  component  flowing  from  high  to  low  bias  terminal  is  much 
larger  than  current  injecting  into  emitter.  The  real  part  of  the  total  small- 
signal  current  remains  approximately  constant  as  can  be  seen  in  figure  4.4(a). 
Since  the  applied  biases  vp (y)  and  Vjy  are  real  in  this  case,  there  is  no  imaginary 
current  flowing  from  high  to  low  bias  terminal.  The  imaginary  currents  at  both 
terminals  decrease  as  vp(y)  decreases  which  results  in  a  decrease  of  the  total 
imaginary  current  (figure  4.4(b)).  Note  the  matched  real  current  and  the 
discrepancy  of  the  imaginary  current  between  the  simulations  and  model.  It 
implies  the  resistance  is  well  modeled  and  the  capacitance  is  not  properly 
described  by  the  uniform  transmission  line  model.  The  reason  is  probably  due 
to  the  assumption  of  the  constant  capacitance  C^  in  the  model,  which  might  not 
be  true  in  the  simulated  device  since  the  intrinsic  base  region  is  not  rectangu- 
lar and  the  doping  density  is  not  uniform  along  the  base-emitter  junction. 
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Figure  4.2.  The  amplitude  of  normalized  base  current  densities  in  the  intrinsic 
base  region  at  different  frequencies.  In  all  cases,  the  small-signal 
voltage  biases  vjsj-vp-  3mV. 
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Figure  4.3.  The  amplitude  of  normalized  base  current  densities  in  the  intrinsic 
base  region  at  different  frequencies.  In  all  cases,  vpj  =  3mV  and  vp 
=  2.7mV.  The  arrows  indicate  the  position  of  minimum  voltage.  For 
the  curves  of  1GHz  and  100MHz,  the  small-signal  voltages  decrease 
monotonically  and  there  is  no  minimum  voltage  inside  the  intrinsic 
region. 
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Figure  4.4(a).  The  real  part  of  the  small-signal  base  current  densities  plotted 
as  a  function  of  bias  vp  when  t>jy  =  3mV.  The  oscillation  frequency 
is  100MHz.  Lines  are  the  results  of  the  model  and  symbols 
indicate  the  data  from  FLOODS  simulations. 
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Figure  4.4(b).  The  imaginary  part  of  the  small-signal  base  current  densities 
plotted  as  a  function  of  bias  vp  when  v^j  -  3mV.  The  oscillation 
frequency  is  100MHz.  Lines  are  the  results  of  the  model  and 
symbols  indicate  the  data  from  FLOODS  simulations. 
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Figure  4.5  illustrates  the  small-signal  j-v  relation  when  the  device  is 
operated  at  10GHz.  The  imaginary  current  shows  the  same  behavior  as  that  at 
the  lower  frequency  operations.  Due  to  the  significant  decrease  of  the  injection 
current  through  the  junction,  the  reduced  amount  of  real  Jp  is  larger  than  the 
increase  of  Jjy;  Thus,  both  the  total  real  and  imaginary  base  currents  decrease 
as  predicted  by  the  simulations  as  well  as  by  model.  Though  there  is  a  numer- 
ical discrepancy  between  the  model  and  simulations,  the  j-v  relation  is  linear 
at  all  frequency  ranges  of  interest,  which  verifies  the  validity  of  eq.  (4-13)  and 
(4-14). 

4.3  Three  Dimensional  A.C.  Distributed  Model 
4.3.1  Model  Development 

The  following  assumptions  are  made  in  the  derivation  of  the  distributed 
model  in  the  base  polysilicon  region:  (i)  The  capacitance  associated  with  the 
polysilicon  base  is  neglected  since  it  is  in  parallel  with  Cd  and  is  much  smaller 
than  Cd.  (ii)  Since  we  limit  our  interests  in  the  biases  before  high  current 
effects  occur,  the  d.c.  voltage  bias  is  approximately  constant  along  FPOLY  as 
can  be  seen  in  figure  3.4.  (hi)  The  admittance  Yd  is  constant  everywhere  in  the 
intrinsic  base  region  due  to  the  small  voltage  variation  in  FPOLY  and  the 
intrinsic  base  region,  (iv)  The  small-signal  voltage  along  lineN  in  NPOLY  is 
constant  due  to  the  small  debiasing  effect  in  NPOLY. 

The  governing  equations  of  the  small-signal  voltage  vPF(y)  and  current 

density  jpp(y)  in  FPOLY  are  given  by 

w  _ 
JPF(y)  =  JPF(W)  +  jJ-^-dy  (4-18) 

y 
and 
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Figure  4.5(a).  The  real  small-signal  base  current  densities  plotted  as  a  function 
of  bias  up  when  i/jy  =  3mV.  The  oscillation  frequency  is  10GHz. 
Lines  are  the  results  of  the  model  and  symbols  indicate  the  data 
from  FLOODS  simulations. 
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Figure  4.5(b).  The  imaginary  small-signal  base  current  densities  plotted  as  a 
function  of  bias  vp  when  vN  -  3mV.  The  oscillation  frequency  is 
10GHz.  Lines  are  the  results  of  the  model  and  symbols  indicate 
the  data  from  FLOODS  simulations. 
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y 

vPF(y)  =  vPF(0)-jjpF(y)pp-dy  (4-19) 

o 
Following  the  same  methodology  used  to  solve  the  d.c.  problem,  the  analytical 

solutions  oi'vPF(y)  and  jpF(y)  can  be  written  as 

vPF(y)  =  ^PF(0)  +  U^^-(ooShaC,(W-y))-cosh(CsW))) 

(4-20) 
jPF(0) 


jpFiy)  ~  sinh(C8W)' 
where 


sinh(Cc(W-y))  (4-21) 


f5P 


ca  =    -¥*■  (4-22) 


'8  -   J  A 


If  the  current  injection  from  the  short  edges  of  emitter  stripe  is  negligible,  then 
the  boundary  conditions  of  small-signal  voltage  uPF(y)  and  current  jPF(y)  are 
given  by 

vPF(0)  =  vN-ApRpeJpF(0)  (4-24) 

where 

C,0n  C,p„Wcoth(C8W) 

C9  =  1  -  -£\  +  C5WRpe  +     5Pp  8  (4-25) 

ApCg  apc8 

Note  that  the  only  variable  in  eq.  (4-23)  and  (4-24)  is  vN  which  is  the  given  bias 

at  lineN. 
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4.3.2  Results  and  Discussion 

The  device  described  in  section  3.2.1  with  an  emitter  stripe  length  of 
lOOjim  is  used  in  the  analysis.  In  all  cases,  d.c.  biases  Vbe=0.83V  and  Vcb=2V 
are  applied  to  the  transistor.  To  improve  the  accuracy  of  the  terminal  j-v 
relations,  the  constants  C5  and  C6  in  equations  (4-13)  and  (4-14)  are  obtained 
directly  from  the  device  simulations.  Figure  4.6  shows  the  amplitude  of  the 
small-signal  voltage  along  FPOLY  at  different  frequencies  while  vN  =  3mV.  At 
low  frequencies,  the  base  current  in  the  intrinsic  base  region  is  dominated  by 
the  current  conducted  through  the  resistance  Rd.  This  current  is  small  and  has 
no  significant  debiasing  effect  in  the  active  and  polysilicon  base  regions.  There 
is  also  no  phase  lag  due  to  the  insignificant  capacitive  characteristics  of  the 
network  as  can  be  seen  in  Fig.  4.7.  The  displacement  current  through  Cd 
increases  following  the  increase  of  the  oscillation  frequency  becoming  more 
important  as  it  dominates  the  base  current  at  higher  frequencies.  The  increas- 
ing current  causes  larger  ohmic  voltage  drops  along  the  FPOLY.  For  instance, 
the  amplitude  of  the  total  small-signal  current  flowing  through  the  FPOLY  at 
10GHz  is  20  times  larger  than  that  at  100MHz.  The  variation  of  small-signal 
voltage  in  FPOLY  is  more  than  72%  at  a  frequency  of  10GHz  whereas  the  varia- 
tion is  within  5%  when  the  operation  frequency  is  100MHz.  The  phase  lag  also 
increases  at  higher  frequencies  due  to  the  prominent  capacitive  characteristics. 
The  modeled  small-signal  base  current  densities  injected  from  the 
polysilicon  contact  to  the  active  base  region  are  shown  in  Fig.  4.8  as  a  function 
of  position  in  the  base  polysilicon.  Figure  4.8(a)  illustrates  the  current  densities 
when  the  operating  frequency  is  100MHz.  At  this  frequency,  the  current  injec- 
tion across  the  junction  is  approximately  constant  since  the  ohmic  voltage  drop 
is  so  small  along  the  base  polysilicon  as  well  as  in  the  intrinsic  base  region. 
However,  the  small  a.c.  debiasing  in  FPOLY  does  cause  small  but  asymmetric 
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Figure  4.6.  The  amplitude  of  the  small-signal  voltage  along  FPOLY  calculated 
by  the  model  when  vjy  =  3mV.  The  emitter  stripe  length  is  lOOum. 
The  voltage  and  current  are  assumed  to  be  symmetric  to  the  center 
of  the  poly. 
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Figure  4.7.  The  phase  lag  of  the  small-signal  voltage  along  FPOLY  computed 
by  the  model.  All  the  phase  lag  refer  to  the  voltage  vN. 
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Figure  4.8(a).  The  small-signal  base  current  densities  injected  from  poly  to 
active  base  region  at  a  frequency  of  100MHz. 
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Figure  4.8(b).  The  small-signal  base  current  densities  injected  from  poly  to 
active  base  region  at  a  frequency  of  10GHz.  Both  the  real  and 
imaginary  total  current  densities  decrease  as  Y  increases. 
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a.c.  biases  between  the  base  contacts.  A  similar  current  dividing  phenomenon 
takes  place  in  the  intrinsic  base  region  as  that  observed  in  the  d.c.  asymmetric 
bias  conditions.  An  obvious  nonuniform  current  distribution  is  observed  in  both 
sides  of  the  base  polysilicon  contacts  and  both  the  real  and  imaginary  currents 
exhibit  the  compensation  property  as  seen  in  the  d.c.  low  current  injection 
cases,  i.e.,  the  current  lost  in  the  FPOLY  is  balanced  by  the  increase  of  current 
in  the  NPOLY,  and  the  total  current  is  approximately  constant  compared  with 
the  case  of  uniform  bias  in  the  FPOLY.  As  the  frequency  increases,  the  debias- 
ing  effect  along  the  FPOLY  becomes  significant  due  to  the  increasing  current 
resulting  in  a  large  decrease  of  a.c.  voltage  across  the  junction  in  the  active 
base  region.  Thus,  the  total  current  decreases  compared  with  the  case  of 
uniform  bias  in  the  FPOLY.  Figure  4.8(b)  shows  the  current  densities  at 
10GHz.  Note  the  a.c.  current  at  10GHz  is  higher  than  that  at  lOOMz  which 
results  in  the  considerable  debiasing  effect  in  FPOLY.  The  decrease  of  real 
current  in  the  FPOLY  is  more  than  the  increasing  current  in  the  NPOLY  and 
the  imaginary  currents  at  both  sides  decrease  toward  the  center  of  polysilicon 
base.  So,  both  the  real  and  imaginary  base  current  calculated  by  a  model  that 
does  not  take  into  account  the  distributed  effects  in  the  polysilicon  base  region 
will  be  overestimated.  Figure  4.8  also  reveals  the  current  redistribution 
modulated  by  the  operating  frequency.  The  ratio  of  current  contributed  by  the 
low  resistive  NPOLY  increases,  and  that  from  the  high  resistive  FPOLY 
decreases  as  the  frequency  increases.  Thus,  in  addition  to  the  intrinsic  base 
region,  the  distributed  effect  in  the  base  polysilicon  region  also  causes  the 
decrease  of  complex  input  impedance  as  the  frequency  increases. 
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Figure  4.9.  The  definition  of  effective  length  Leg-used  in  this  work.  The  shaded 
areas  indicate  the  total  current  1^  f  associated  with  each  poly  base, 
which  is  also  equivalent  to  the  product  of  Lgffj^p  and  JjgF(0). 
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4.3.3  Application  and  Verification 

To  characterize  the  high  frequency  performance  of  a  transistor,  2D 
device  simulations  are  frequently  used.  A  key  assumption  of  the  simulations  is 
that  the  voltage  and  current  distribution  in  the  third  dimension  are  uniform. 
The  results  obtained  using  our  3D  distributed  model  suggest  the  assumption 
may  be  invalid  in  certain  structures,  for  some  cases,  a  3D  device  simulator  is 
required.  However,  3D  device  simulations  are  rarely  used  in  device  design  due 
to  extensive  demands  on  hardware  and  CPU  time.  The  use  of  a  2D  device 
simulator  coupled  with  the  analytical  model  to  handle  the  variation  of  voltage 
and  current  in  the  third  dimension  is  an  efficient  way  to  characterize  these 
transistors. 

The  distributed  model  is  coupled  to  FLOODS  by  effective  contact 
lengths.  Referring  to  the  base  current  density  plot  as  shown  in  Fig.  4.8,  the 
total  current  I^p  associated  with  each  side  of  the  polysilicon  base  contact  is  the 
integration  of  the  current  density  J^pfy)  along  the  polysilicon  length  (y-direc- 
tion)  which  can  be  equivalent  to  the  current  density  at  the  edge  of  polysilicon 
contact,  J^p(0),  times  the  effective  contact  length  Leff^p,  i.e. 
IN  F  -  JN  F(0)  xLep  F.  Figure  4.9  illustrates  the  definition  of  effective 
lengths  in  the  NPOLY  and  FPOLY  for  the  current  densities  shown  in  Fig. 
4.8(a).  Leff  of  NPOLY  is  longer  than  the  actual  polysilicon  contact  length  Lactuai 
since  the  current  density  increases  along  NPOLY.  In  contrast,  Leff  of  FPOLY  is 
shorter  than  Lactuai  due  to  the  decreasing  current  density  along  FPOLY.  Thus, 
if  we  obtain  the  edge  current  densities  J^piO)  from  FLOODS  simulations  and 
multiply  by  the  Lefwp  from  the  model,  we  get  the  actual  contact  current  from 
the  nonuniform  3D  effect  of  the  transistor. 

Figure  4.10  illustrates  the  implementation  of  a  realistic  3D  structure,  as 
shown  in  figure  2.14  and  2.15,  to  a  simulated  structure.  The  differences 
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Figure  4.10.  Schematic  diagram  of  a  BJT  with  parasitics  used  in  the  asymmetric 
simulations.  The  total  current  flowing  through  each  contact  of  the  active 
device  IN  F  equals  the  current  density  JNF  times  the  effective  length  LeffN  F 
of  the  contact.  The  connection  resistance  Rpe  between  NPOLY  and  FPOLY 
is  6.3Q. 
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Figure  4.11.  The  comparison  of  measured  and  simulated  Sn  as  a  function  of 
frequency.  The  asymmetric  simulation  uses  the  structure  shown 
in  figure  4.10  and  the  symmetric  simulation  utilizes  the  structure 
shown  in  figure  2.16.  The  cutoff  frequency  fp  of  this  device  is  10.3 
GHz. 
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between  figure  4.10  and  2.16  are  that  (i)  a  lumped  resistor  Rpe  is  placed 
between  the  NPOLY  and  FPOLY  contacts  to  depict  the  connection  resistance 
between  them,  (ii)  Leff  's  are  used  in  the  base  contacts  to  describe  the  nonuni- 
form effect  in  each  contact.  Due  to  the  ohmic  voltage  drop  across  Rpe,  the  bias 
at  FPOLY  is  lower  than  the  bias  at  NPOLY  which  properly  represents  the 
actual  asymmetric  bias  condition  at  the  edge  of  the  poly  bases. 

To  characterize  the  high  frequency  performance  of  a  transistor,  scatter- 
ing parameters  (s-parameters)  are  often  used  because  they  are  easier  to 
measure  at  high  frequencies  than  other  parameters  [Sze81].  However,  it  is 
difficult  to  properly  describe  the  transmission  line  with  finite  RLC  circuit 
elements  in  the  simulations;  therefore,  y-parameters  are  extracted  from  device/ 
circuit  simulations  and  then  converted  to  equivalent  s-parameters  [Gon84]. 
Here  we  will  only  focus  on  the  input  reflection  coefficient  Sn  or  input  imped- 
ance Zin  of  the  device  since  the  base  distributed  effects  have  a  direct  relation 
with  Sn  and  Zin  in  a  common-emitter  configuration.  Sn  and  Zin  are  related  by 

a  Zin-Zol  (4.26) 

117,7 

where  Zin  is  the  input  impedance  of  the  transistor  which  is  associated  with  the 
base  and  emitter  characteristics  in  a  common-emitter  configuration,  and  Zol  is 
the  complex  characteristic  impedance  of  the  transmission  line  which  is  given 
by  [Gon84] 


R0  +  jwLc 


Z     =         "      ,       "  (4-27) 

01       ^Go  +  jwC0 

where  R0,  G0,  L0,  and  C0  are  the  resistance,  conductance,  inductance  and 
capacitance  per  unit  length  of  the  transmission  line.  At  microwave  frequencies, 
R0  and  G0  are  usually  negligible  and  the  transmission  line  is  said  to  be  lossless. 

!l~0 

In  a  lossless  transmission  line,  Zol  =    j-^-  and  is  a  real.  If  assuming  the  induc- 
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tive  component  of  a  transistor  is  insignificant  at  the  frequency  range  of 
interest,  then  Zin  can  be  represented  by  a  series  RC  equivalent  circuit.  At  a 

specific  angular  frequency  w,  Z-     =  9?(u>)  +  -: — -  .  Substituting  into  eq.  (4- 

26),  it  gives 

(iR{w)+zoo-(R(w)-zo0+-1^yj-^) 

Sn  =  WUJ£) (4.28) 

{R{w)  +  Zol)2  +  — - 

w  C(w) 


The  comparison  of  the  measured  and  simulated  Sns  are  shown  in  Fig. 
4.11  as  a  function  of  frequency.  The  asymmetric  simulation  uses  the  structure 
shown  in  Fig.  4.10  with  the  effective  contact  lengths  in  the  device/circuit  inter- 
face which  take  into  account  the  nonuniform  effects  in  the  base  polysilicon 
region.  The  symmetric  simulation  utilizes  the  same  structure  with  the  actual 
contact  lengths  in  the  interface  and  Rpe=Q  which  results  in  a  symmetric  bias  at 
the  NPOLY  and  FPOLY.  For  the  intrinsic  simulation,  the  circuit  network  is  not 
included.  It  is  obvious  that  the  simulation  cannot  predict  correctly  if  the 
parasitics  are  not  being  accounted  for.  The  real  part  of  Sn  is  reasonably 
described  by  both  the  asymmetric  and  symmetric  simulations  at  the  low 
frequencies.  As  the  frequency  increases,  the  result  of  the  symmetric  simulation 
deviates  from  the  measured  data,  whereas  the  asymmetric  simulation  result 
remains  in  agreement  with  the  measured  one.  For  imaginary  Sn,  there  is  a 
discrepancy  between  measurement  and  symmetric  simulation,  while  the 
asymmetric  simulation  result  agrees  with  the  measured  data  over  the  low 
frequency  range,  and  as  the  frequency  increases  the  measured  results  deviate 
from  the  asymmetric  simulations. 

An  analysis  of  the  input  impedance  Zin  reveals  that  91  (w)  is  a  weak 
function  of  frequency,  but  %  (w)  is  a  strong  function  of  frequency.  From  Eq.  (4- 
28),  it  is  clear  to  see  that  the  w%  (w)  term  dominates  the  real  part  of  Sn  at  low 
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frequencies  and  the  Si(w)  term  becomes  as  important  as  w%(w)  at  intermedi- 
ate frequencies.  At  high  frequencies,  the  real  Sjj  is  dominated  by  the  SR(w) 
term.  The  imaginary  part  of  Sn  is  governed  by  the  w%(w)  term  over  all 
frequency  ranges.  From  these  results  we  can  conclude  that  the  resistive  charac- 
teristics of  the  uniform  transmission  line  model  in  the  active  base  region  is 
properly  modeled  if  the  distributed  effect  in  the  polysilicon  base  region  is 
included.  The  model  can  also  accurately  describe  the  capacitive  property  for 
low  and  medium  frequency  applications  and  becomes  less  accurate  at  high 
frequencies.  This  discrepancy  may  be  due  to  the  nonuniform  capacitance  in  the 
intrinsic  base  region  which  in  fact  is  not  rectangular  resulting  in  a  nonuniform 
doping  profile  in  the  base-emitter  junction  along  the  x-direction. 

I  Sn  |  2  represents  the  ratio  of  the  power  reflected  from  input  port  to  the 
power  available  at  the  input  port  and  is  associated  with  the  maximum  avail- 
able power  gain  and  stability  factor  of  a  transistor.  At  the  cutoff  frequency  fy  (= 
10.3  GHz),  the  |  Sn  |  2  predicted  by  the  asymmetric  and  symmetric  simulations 
is  higher  than  the  measured  data  by  3%  and  32%,  respectively. 

4.4  Summary 
An  a.c.  model  has  been  developed  to  study  the  high  frequency  distrib- 
uted effect  of  voltage  and  current  in  the  intrinsic  and  polysilicon  base  regions. 
The  small-signal  voltage  and  current  distributions  in  both  regions  have  been 
described  completely  by  the  model.  By  the  introduction  of  the  effective  contact 
length,  the  model  has  been  applied  to  a  2D  device  simulator  for  quasi-3D 
simulations.  Sn  was  extracted  from  different  simulations  and  compared  with 
the  measured  data.  The  agreement  between  measured  and  simulated  S1:L  has 
been  improved  by  the  consideration  of  the  nonuniform  effect  in  the  third 
dimension  of  the  device.  A  modification  of  the  nonuniform  capacitive  character- 
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istic  in  the  small-signal  network  is  required  to  further  improve  the  accuracy  of 
the  imaginary  S^  at  high  frequencies. 


CHAPTER  5 
VERIFICATION  OF  DC  DISTRIBUTED  MODEL  AND 
GEOMETRY  DEPENDENT  DEVICE  PERFORMANCE 
DUE  TO  DISTRIBUTED  EFFECTS 


5.1  Introduction 


Stripe  and  interdigitated  configurations  (Fig.  5.1)  are  two  popular  layout 
geometries  used  in  high  speed  BJT  design.  Designers  usually  assume  the 
current  and  voltage  distribution  in  the  z-direction  of  the  polysilicon  contacts 
are  uniform,  and  the  large-  and  small-signal  currents  are  linearly  proportional 
to  the  length  of  the  stripes  or  fingers.  The  analysis  of  Ch.  3  and  Ch.4  shows  this 
assumption  is  not  correct  for  striped  single  metal  base  contact  structure  (SMB) 
in  some  applications.  Similar  debiasing  and  distributed  effects  also  take  place 
in  the  fingers  of  interdigitated  structures.  However,  DC  behaviors  of  varying 
types  of  layouts  and  sizes  can  differ  due  to  varying  degrees  of  influence  by 
distributed  effects  on  the  different  structures. 

The  influence  of  the  distributed  effects  on  the  structures  strongly 
depends  on  the  sheet  resistance  of  the  base  polysilicon  contacts.  A  current  path 
of  low  resistance  reduces  debiasing  and  diminishes  distributed  effects. 
Silicided  polysilicon  technologies  are  often  utilized  in  current  high  speed 
semiconductor  devices  to  reduce  the  resistivity  and  RC  delay.  The  sheet  resis- 
tance of  contemporary  silicide  technologies  ranges  from  single  to  teens  Q/sq 
[W0I86,  Chi91,  Yam93].  A  low  resistance  silicide  process  usually  increases  the 
complexity  of  fabrication  and  cost.  Thus,  the  adequate  choice  of  the  sheet  resis- 
tance of  polysilicon  contacts  is  one  of  the  key  issues  in  the  device  design. 
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Figure  5.1  Three  layout  structures  for  high  speed  BJTs.  (a)  Double  Metal  Base 
contact  structure  (DMB).  (b)  Single  Metal  Base  contact  structure 
(SMB).  (c)  Interdigitated  structure. 
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Three-dimensional  simulation  capability  has  been  implemented  into 
FLOODS  allowing  a  direct  study  of  the  3D  distributed  effects.  In  this  chapter, 
the  3D  simulations  are  demonstrated  and  some  questions  about  3D  simula- 
tions are  discussed.  The  results  of  the  DC  distributed  model  derived  in  Ch.  3 
are  compared  with  results  from  3D  simulation  to  verify  the  accuracy  and 
limitations  of  the  model.  The  device  performance  of  double  polysilicon  BJTs 
with  DMB,  SMB,  and  interdigitated  layout  structures  is  studied  for  different 
geometry  sizes  and  polysilicon  contact  sheet  resistances. 

5.2  3D  Device  Simulations  and  Simulation  Structures 

The  box's  (finite  difference)  method  is  used  as  the  discretization  scheme 
for  the  3D  simulations.  The  algorithm  to  solve  a  3D  problem  is  the  analogous 
as  that  of  2D  simulations.  Using  the  discretized  form  of  Poisson  equation  as  an 
example, 

^Dij-dSij*  prdVt  (5-1) 

j 

where  Dij  is  the  potential  flux  along  edge  ij  between  nodes  i  andj,  and  p;  is  the 

net  charge  associated  with  node  i.  In  3D  simulations,  dSy  is  the  cross  section 

area  bisecting  edge  ij,  and  dV^  is  the  volume  associated  with  node  i,  while  dStj 

is  the  cross  section  length  bisecting  edge  ij,  and  dV^  is  the  area  related  to  node 

i  in  2D  simulations.  For  object-oriented  languages  like  C++,  the  program  can 

automatically  choose  the  proper  dimensions  for  dSy  and  dVi  in  2D  or  3D 

simulations.  Therefore,  only  a  small  amount  of  additional  coding  is  necessary 

to  accommodate  the  3D  capability  into  FLOODS. 

One  main  difference  between  the  2D  and  3D  simulations  is  the  density 

of  the  sparse  matrix.  For  a  node  in  a  ID  mesh,  there  are  two  adjacent  nodes. 
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These  adjacent  nodes  become  six  in  a  2D  right  triangular  mesh  and  increase  to 
14  in  a  3D  tetrahedral  mesh.  Each  adjacent  node  represents  a  nonzero  entry  in 
one  row  of  the  matrix.  Therefore,  the  sparse  matrix  of  3D  simulations  is  denser 
than  the  matrix  of  2D  simulation  for  the  same  number  of  nodes.  For  example, 
the  number  of  nonzero  entries  of  a  drift-diffusion  2D  simulation  with  2160 
nodes  is  99218.  For  a  similar  3D  simulation  with  2080  nodes,  the  number  of 
nonzero  entries  becomes  163294  which  requires  twice  as  much  memory  and 
CPU  time  as  compared  with  2D  simulations. 

The  3D  simulation  structure  of  a  high  speed  npn  BJT  is  illustrated  in 
Fig.  5.2.  The  material  of  the  bottom  cube  is  silicon  and  the  two  sections  stacked 
on  the  silicon  are  p-type  base  polysilicon  contacts.  For  a  DMB  layout  structure, 
the  metal  contacts  are  located  at  faces  A  and  B  and  are  connected  to  a  voltage 
source.  The  metal  contacts  of  a  SMB  structure  are  made  on  face  A  and  D  where 
the  left  section  of  polysilicon  acts  as  NPOLY  and  the  right  section  of  polysilicon 
is  FPOLY.  There  is  a  resistor  connected  between  the  two  contacts  to  represent 
the  resistance  of  polysilicon,  Rpe,  between  NPOLY  and  FPOLY,  and  Contact  A 
is  connected  to  a  voltage  source.  The  metal  contacts  of  the  interdigitated  struc- 
ture are  made  on  faces  C  and  D  and  are  connected  to  a  voltage  source.  The 
sheet  resistance  of  the  polysilicon  contact  R^  is  controlled  by  the  specified 
constant  hole  mobility  in  the  polysilicon  material  which  is  independent  of  the 
mobility  in  the  silicon  region. 

The  accuracy  and  convergence  capabilities  of  device  simulations 
strongly  depend  on  the  mesh  density  and  location  associated  with  the  device 
structure.  The  regions  where  large  variation  of  doping,  potential  and  carrier 
densities,  etc.,  exist  require  fine  grid  to  properly  describe  these  physical  param- 
eters at  the  specific  positions  of  the  regions.  Figure  5.3  illustrates  the  Gummel 
plot  of  a  5um  stripe  length  npn  BJT  for  varying  mesh  densities  in  the  3D 
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Figure  5.2  The  3D  simulation  structure  of  a  BJT.  A,  B,  C  and  D  mark  the  faces 
where  a  metal  contact  is  attached  to  the  polysilicon. 
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Figure  5.3  The  Gummel  plots  of  a  BJT  simulated  with  a  varying  number  of  grid 
points  in  the  mesh  structure. 
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simulations.  The  first  mesh  consists  of  about  400  nodes  in  each  xy-plane  with 
4  xy-planes  (3  segments)  distributed  at  5um  intervals  along  the  z-direction. 
There  are  1666  nodes  in  this  mesh.  The  second  mesh  refines  the  grids  in  the  x- 
direction  and  the  number  of  nodes  becomes  2032.  Due  to  the  finer  grid  in  the 
base-collector  depletion  region,  the  second  mesh  predicts  a  wider  base  width. 
Thus,  the  collector  current  of  the  finer  mesh  structure  is  9%  smaller  at 
Vbe=0.7V  compared  with  the  1666  node  mesh.  The  base  currents  of  these  two 
meshes  do  not  show  significant  differences  since  the  doping  density  and  deple- 
tion region  boundaries  around  the  base-emitter  junction  is  not  significantly 
changed  by  the  finer  mesh.  The  third  mesh  further  scales  the  grid  in  the  y- 
direction.  In  this  case,  the  base  doping  at  the  sidewall  of  the  base-emitter 
junction  decreases,  reducing  the  injection  of  holes  into  the  emitter.  Therefore, 
the  base  current  is  decreased  by  25%  and  20%  at  V^e=0.TV  compared  with  the 
first  and  second  meshes,  respectively.  Since  few  electrons  are  injected  from  the 
highly  doped  sidewall  into  the  base,  the  collector  current  remains  unchanged. 
Further  scaling  of  the  mesh  in  the  xy-planes  improves  the  accuracy,  but 
this  improvement  becomes  negligible  after  the  mesh  is  fine  enough.  If  the 
doping  profile  in  the  z-direction  is  uniform,  the  accuracy  of  the  device  simula- 
tions is  insensitive  to  the  mesh  density.  For  example,  the  use  of  three  segments 
or  five  segments  of  mesh  does  not  show  significant  difference  for  the  I-V  curve 
of  a  SMB  device  where  a  significant  debiasing  occurs  in  the  z-direction.  This  is 
probably  due  to  the  smooth  change  of  voltage  which  can  be  properly  approxi- 
mated by  these  segments.  However,  it  is  observed  that  the  use  of  a  very  large 
grid  size  in  the  z-direction  reduces  the  convergence  capability  due  to  the  large 
aspect  ratio  in  the  tetrahedral  unit  of  the  mesh.  Therefore,  if  the  piece-wise 
linear  voltage  can  properly  approximate  the  real  voltage  in  the  z-direction  and 
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convergence  is  not  a  problem,  the  grid  size  in  the  z-direction  can  be  set  as  large 
as  possible. 

Limited  by  the  available  computer  resources  (Sun  Ultrasparc  2  worksta- 
tions with  1.6G  accessible  memory),  all  the  3D  simulations  in  this  and  the 
following  chapters  use  the  first  mesh  (400  nodes  in  xy-plane  as  shown  in  Fig. 
5.2).  The  grid  size  in  the  z-direction  varies  from  one  to  2.5|im  depending  on  the 
length  of  the  device.  The  total  number  of  nodes  in  the  meshes  is  limited  to  2100 
nodes  which  requires  about  1.5G  of  memory  for  AC  simulations.  Though  the 
loose  mesh  in  the  xy-plane  introduces  error  in  the  I-V  curves,  the  trend  and  the 
characteristics  of  the  distributed  effects  in  the  transistors  can  still  be  properly 
described  by  this  3D  simulations. 

Usually  a  base  polysilicon  contact  covers  over  both  the  base  diffusion 
region  and  the  field  oxide  region.  The  sheet  resistance  on  the  top  of  field  oxide 
is  typically  higher  than  that  on  the  diffusion  region.  To  simplify  the  simulated 
structure,  the  nonuniform  sheet  resistance  in  the  y-direction  can  be  repre- 
sented by  a  single  effective  sheet  resistance  Rs^  which  is  obtained  by 
simulating  a  complete  3D  polysilicon  structure  including  the  base  diffusion  and 

field  oxide  region:  the  resistance  of  polysilicon  R  is  computed  from  the  terminal 

W  ol 
voltage  and  current,  then  Rsh  =  R  ■      po  - ,  where  Wpoiy  and  Lpoiy  are  the 

poly 

width  and  length  of  the  polysilicon,  respectively.  This  effective  sheet  resistance 
can  properly  describe  the  debiasing  effects,  although  the  simulation  structure 
as  shown  in  Fig.  5.2  does  not  include  the  entire  polysilicon  structure. 

5.3  DC  Characteristics 

As  discussed  in  section  3.3.3,  FPOLY  in  the  SMB  structure  represents  a 
highly  resistive  current  path  and  thus  has  a  bigger  effective  base  resistance 
than  the  DMB  structure.  For  the  same  reason,  both  fingers  of  the  interdigitated 
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Figure  5.4.  (a)  The  Gummel  plots  of  BJTs  with  DMB  and  SMB  layout 
structures  obtained  by  3D  simulations.  Four  differentials  of  base 
polysilicon  are  used  in  the  simulations  for  the  SMB  structure. 
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Figure  5.4.  (b)  The  simulated  Gummel  plots  of  BJTs  with  different  layout 
structures.  The  value  of  Rs^  in  polysilicon  contacts  is  12.5H/sq. 
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structure  have  highly  resistive  current  paths  compared  with  the  paths  in  the 
DMB  device.  Therefore,  the  DMB  structure  has  the  lowest  effective  base  resis- 
tance, and  the  interdigitated  structure  has  the  highest  one.  The  effective  base 
resistance  can  be  calculated  from  the  deviation  of  base  current  from  the  ideal 
exponential  curve  in  a  Gummel  plot.  The  difference  of  the  voltage  is  mainly  due 
to  the  ohmic  voltage  drop  in  the  base  and  emitter  regions  as  expressed  in  Eq. 
(2-18).  Qualitatively,  a  larger  degree  of  base-current-curve-bending  in  a 
Gummel  plot  indicates  a  larger  base  resistance. 

The  Gummel  plots  of  a  5fim  striped  B  JT  with  different  layout  structures 
and  i2s/jS  are  shown  in  Fig.  5.4.  At  low  biases  (<  0.85V),  the  debiasing  is  very 
small  either  in  the  active  base  or  in  polysilicon  base  contacts.  The  base  current 
increases  exponentially  with  V^  and  is  independent  of  layout  structures  and 
Rsh.  As  the  bias  increases,  the  base  current  curves  for  different  layouts  andi?s/j 
depart  from  each  other.  It  is  clear  that  the  base  resistance  of  the  SMB  structure 
is  always  larger  than  that  of  the  DMB,  and  a  larger  sheet  resistance  increases 
the  difference.  For  the  same  value  of  Rsh,  the  degree  of  bending  for  the  inter- 
digitated structure  is  larger  than  that  of  the  SMB  transistor  as  plotted  in  Fig. 
5.4(b).  Therefore,  the  layout  with  the  DMB  structure  achieves  the  lowest  base 
resistance,  and  the  interdigitated  structure  gives  the  highest  base  resistance. 

Figure  5.5  shows  the  hole  Fermi  potential  along  FPOLY  of  a  SMB  device 
at  V"5e=0.95V  computed  by  3D  simulations  and  models.  The  current  density  in 
the  xy-plane  of  FPOLY  in  the  model  is  assumed  to  be  constant.  However,  the 
current  density  in  FPOLY  of  the  3D  simulations  is  not  uniform  due  to  the 
nonuniform  doping  density  in  xy-plane.  Thus,  for  the  purpose  of  comparison, 
two  potential  curves  are  extracted  from  the  top  surface  of  FPOLY,  one  from  the 
inner  side  (y=0.45|im)  and  the  other  from  the  outer  side  (y=0.65|im).  As  seen  in 
Fig.  5.5,  the  potential  in  FPOLY  descends  monotonically  from  the  contact  side 
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Figure  5.5  The  potential  in  FPOLY  of  a  SMB  structure  at  a  bias  of  0.95V.  The 
metal  contact  is  located  at  z=0  which  is  attached  to  a  voltage 
source  by  a  resistor  Rpe.  The  inner  potential  is  extracted  from  the 
inner  side  (y=0.45um)  of  the  polysilicon  top  surface  in  the  3D 
simulation  structure.  The  outer  potential  is  from  the  outer  side  of 
the  polysilicon  top  at  y=0.65(im. 
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(z=0)  toward  the  center  (z=5|im),  and  a  higher  R^  causes  a  larger  debiasing. 
There  is  also  a  connecting  resistance  Rpe  between  base  polysilicon  contacts 
resulting  in  additional  ohmic  voltage  drop  at  the  edge  of  FPOLY.  For  instance, 
the  potential  difference  between  the  edge  and  center  of  a  FPOLY  with  #^=12.5 
QJsq  is  about  lOmV  and  voltage  drop  across  Rpe  is  about  4mV  obtained  from  the 
simulations.  This  debiasing  results  in  a  significant  decrease  of  the  base  current 
compared  with  the  DMB  structure  as  seen  in  Fig.  5.4.  The  curves  predicted  by 
the  model  follow  the  same  trend  and  basically  fall  within  the  range  of  the 
simulations.  Though  there  is  a  discrepancy,  the  current  from  each  contact 
computed  by  the  model  is  about  the  same  as  that  obtained  by  the  3D 
simulations. 

From  the  analysis  in  section  3.2.4,  the  biases  in  two  base  polysilicon 
contacts  of  the  DMB  structure  are  constant  along  z-direction  and  equal  to  each 
other.  The  debiasing  simply  takes  place  in  y-direction  symmetrically  and  the 
currents  flowing  into  each  base  polysilicon  contact  are  the  same.  For  a  SMB 
structure,  the  voltage  in  NPOLY  is  uniform  but  voltage  decreases  along  FPOLY 
resulting  in  asymmetric  biases  in  the  two  polysilicon  base  contacts.  Conse- 
quently, the  currents  from  two  contacts  become  asymmetric  as  well.  Figure  5.6 
illustrates  the  simulated  potential  profile  in  the  intrinsic  base  region  beneath 
the  emitter  of  a  SMB  structure  with  Rs^=12.5  QJsq.  It  is  observed  that  debias- 
ing in  the  intrinsic  base  region  results  in  a  voltage  decreasing  toward  the 
center  in  the  y-direction.  The  voltage  decrease  in  FPOLY  also  causes  the  poten- 
tial to  decrease  in  the  z-direction.  Consequently,  the  minimum  voltage  point 
moves  to  the  right  (positive  y)  and  toward  the  center  of  the  stripe  (Z=5um)  due 
to  the  potential  decrease  in  FPOLY.  Figure  5.6  shows  exactly  the  same  charac- 
teristics as  Fig. 3. 5  computed  by  the  model. 
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Figure  5.6  The  potential  in  the  intrinsic  base  region  of  a  SMB  transistor 
obtained  from  the  3D  simulations  at  a  bias  of"Vje=l,lV. 
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Figure  5.7  A  comparison  between  simulated  and  modeled  base  current  ratio  as 
a  function  of  Vf,e  at  four  different  values  of  J2S^. 
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The  asymmetric  biases  at  the  base  polysilicon  contacts  cause  the 
movement  of  the  minimum  voltage  point  in  the  intrinsic  base  region  resulting 
in  unequal  current  contributions  from  NPOLY  and  FPOLY.  The  ratio  y  of  the 
current  contributed  by  NPOLY  vs.  current  from  FPOLY  of  SMB  structures  is 
greater  than  unity  since  the  voltage  at  NPOLY  is  always  higher  than  the 
voltage  of  FPOLY.  Fig.  5.7  illustrates  the  simulated  and  modeled  y  as  a  function 
of  bias  and  Rg^.  For  a  very  low  R^  device  (<  1Q),  the  ratio  y  approaches  unity 
since  the  voltage  difference  between  these  two  polysilicon  contacts  is  negligible 
even  at  high  biases.  The  debiasing  is  also  not  significant  in  the  devices  with 
intermediate  Rs^  (1Q  <  R^  <  150)  at  a  bias  less  than  0.9V.  The  further  increase 
of  V^g  causes  larger  ohmic  voltage  drops  along  FPOLY  and  enlarges  the  voltage 
difference  between  NPOLY  and  FPOLY  which  enhances  the  movement  of  the 
minimum  voltage  point  toward  FPOLY  side  in  the  intrinsic  base  region.  Conse- 
quently, the  current  from  NPOLY  increases,  and  the  current  from  FPOLY 
decreases,  while  the  ratio  y  becomes  much  greater  than  one.  The  increase  ofi?s^ 
tends  to  trigger  the  movement  of  the  minimum  voltage  point  at  a  lower  bias 
and  induces  larger  movement  at  the  same  bias  compared  with  the  lower  R^ 
devices.  Thus,  a  very  unbalanced  current  contribution  from  the  two  polysilicon 
contacts  is  observed.  Since  the  current  contribution  from  FPOLY  is  negligible 
at  high  bias  for  larger  values  of  Rsh,  the  performance  of  the  SMB  device  is 
similar  to  that  of  a  single  base  polysilicon  contact  device. 

Figure  5.7  shows  that  the  model  predicts  the  current  ratio  y  accurately 
for  devices  with  low  and  intermediate  R^.  The  modeled  y  is  also  well  matched 
to  the  simulation  results  of  devices  with  higher  Rsh  until  it  goes  into  the  high 
current  region.  In  this  region,  the  model  overestimates  the  ratio  y  The  error 
comes  from  the  failure  of  the  linearized  assumption  in  Eq.  (3-10)  where  the 
voltage  difference  between  two  base  polysilicon  contacts  is  assumed  to  be  much 
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Figure  5.8  The  base  current  plotted  as  a  function  of  base  polysilicon  width  and 
sheet  resistance  obtained  from  3D  simulations. 
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smaller  than  two  times  the  thermal  voltage  V0  (~52mV  at  300K).  For  those 
devices  with  large  Rs^  operating  in  the  high  current  regime,  the  voltage  differ- 
ence between  the  base  polysilicon  contacts  approaches  2V0.  Under  this 
condition,  the  increase  of  JjvPOLY  following  the  decreasing  V^is  less  than  that 
predicted  by  the  linearized  model  due  to  the  additional  voltage  drop  in  the  base 
region.  Consequently,  there  is  less  IjqpoLY  in  tne  rea^  situation  resulting  in  the 
overestimation  of  yby  the  model. 

Figure  5.8  shows  the  base  currents  of  SMB  and  DMB  structures  with 
different  polysilicon  contact  widths  Wpoiy  The  wider  contact  reduces  the  resis- 
tance and  decreases  the  ohmic  voltage  drop  along  the  base  current  path.  Thus, 
the  base  current  (and  also  the  collector  current)  of  devices  with  wider  Wpoiy  is 
slightly  larger  than  that  of  a  narrower  Wpoiy  However,  the  use  of  wider  Wpoiy 
does  not  significantly  eliminate  the  distributed  effects.  Rsh  is  still  the  dominant 
parameter  in  determining  the  debiasing  and  distributed  effects.  Though  a 
narrow  Wpoiy  is  used  for  the  devices  studied  in  Ch.  5  and  6,  the  same  results 
are  applicable  to  wider  Wpoiy  devices. 

A  circuit  designer  usually  assumes  the  current  driving  capability 
(current  per  unit  length)  is  linearly  proportional  to  the  length  of  the  device 
stripe.  If  twice  the  amount  of  current  driving  is  required,  the  layout  stripe 
length  is  simply  doubled.  Figure  5.9  plots  the  base  current  driving  capability 
as  a  function  of  stripe  length  for  various  layout  structures  and  sheet  resis- 
tances. The  base  current  of  the  DMB  structure  is  seen  to  be  linearly 
proportional  to  the  stripe  length.  For  the  SMB  and  interdigitated  structures, 
the  unity-length  current  driving  capability  degrades  as  the  stripe  length 
increases,  and  the  larger  Rsh  causes  more  significant  decay  due  to  a  larger 
voltage  decrease  in  the  base  polysilicon.  For  the  same  Rsh,  the  degradation  of 
the  interdigitated  structure  is  larger  than  the  SMB  structure  since  both  base 
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Figure  5.9  The  base  current  vs.  stripe  length,  (a)  Comparison  between  DMB 
and  SMB  structures,  (b)  Comparison  between  DMB  and 
interdigitated  structures. 
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polysilicon  contacts  of  the  interdigitated  structure  suffers  from  the  debiasing 
effects.  The  collector  current  driving  capability  also  shows  the  same  character- 
istics as  that  of  base  current  illustrated  in  Fig.  5.9.  Since  the  input  and  output 
current  driving  is  very  sensitive  to  R^,  the  silicided  process  is  necessary  and 
very  important  to  the  DC  characteristics  of  the  DMB  and  interdigitated 
structures. 

5.4  Summary 

The  base  resistance  of  BJTs  is  shown  to  strongly  depend  on  the  3D 
distributed  effects.  Due  to  different  degrees  of  influence  by  the  effects,  DMB, 
SMB  and  interdigitated  layout  structures  present  different  I-V  relations  and 
have  different  values  of  base  resistance.  The  current  driving  capability  is  also 
affected  by  the  distributed  effects.  The  DMB  layout  is  shown  to  have  the  best 
current  driving  capability  and  lowest  base  resistance,  while  the  interdigitated 
structure  has  the  largest  base  resistance  and  the  worst  driving  capability  for 
the  same  process. 

Three-dimensional  device  simulations  have  been  used  to  verify  the 
accuracy  and  limitation  of  the  quasi-3D  distributed  model.  The  results  show 
the  model  is  valid  if  the  voltage  difference  between  the  two  polysilicon  base 
contacts  is  smaller  than  two  times  the  thermal  voltage.  The  condition  is 
violated  only  if  the  device  has  a  large  polysilicon  sheet  resistance  and  operates 
at  high  biases  where  there  is  less  interest  in  practical  device  applications. 


CHAPTER  6 
VERIFICATION  OF  AC  DISTRIBUTED  MODEL  AND 
GEOMETRY  DEPENDENT  AC  DEVICE  PERFORMANCE 
DUE  TO  DISTRIBUTED  EFFECTS 


6.1  Introduction 

The  small-signal  current  flowing  in  the  intrinsic  and  polysilicon  base 
regions  of  a  BJT  also  represents  the  distributed  characteristics  and  build-up  of 
small-signal  debiasing  along  the  current  paths.  The  input  conductance  of 
common-emitter  BJTs  increases  following  the  increase  of  operating  frequency 
due  to  the  larger  input  small-signal  current  flowing  through  the  capacitance 
associated  with  base.  Therefore,  more  significant  debiasing  effects  are  observed 
at  higher  frequencies.  For  a  SMB  structure,  the  increasing  small-signal  current 
creates  larger  asymmetric  biasing  resulting  in  unbalanced  small-signal 
current  drive  in  the  two  polysilicon  base  contacts.  For  an  interdigitated  struc- 
ture, the  biases  in  the  base  contacts  are  symmetric  and  the  debiasing  takes 
place  in  both  base  fingers.  Thus,  the  use  of  various  layout  geometries  can  result 
in  different  high  frequency  characteristics  for  the  same  device  process. 

In  this  chapter,  the  AC  distributed  model  derived  in  Ch.  4  and  the 
results  of  quasi-3D  simulations  are  verified  by  comparing  them  with  the  3D 
simulation  results.  The  use  of  various  sheet  resistances  R^  in  the  base  polysil- 
icon contacts  on  the  AC  performance  of  BJTs  is  investigated.  The  AC 
characteristics  of  different  layout  structures  and  geometry  sizes  are  also 
studied  based  on  the  3D  simulations  and  the  importance  of  the  distributed 
effects  on  the  RF  circuit  design  is  demonstrated  as  well. 
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6.2  Verification  of  AC  Model  and  Quasi-3D  Simulations 

The  AC  performance  of  B JTs  with  various  layout  structures  are  studied 
in  this  chapter  and  characterized  at  biases  of  Vbe=0.85V  and  Vcb=2V  where 
the  base-push  out  effect  has  not  taken  place,  and  DC  debiasing  does  not  signif- 
icantly affect  the  base  and  collector  currents  (refer  to  Fig.  5.4).  Figure  6.1  shows 
the  modeled  small-signal  voltage  in  FPOLY  of  a  SMB  structure  characterized 
at  four  frequencies  when  a  small-signal  voltage  is  applied  to  the  base  voltage 
source.  It  is  observed  that  the  magnitude  of  the  small-signal  voltage  decays 
along  FPOLY  and  the  debiasing  is  larger  at  higher  frequencies  due  to  the 
increasing  small-signal  current  level.  There  is  also  a  phase  lag  in  the  voltage 
along  FPOLY  and  the  degrees  of  delay  following  the  increase  in  frequency.  At 
very  high  frequencies,  the  phase  delay  near  the  front  side  of  FPOLY  is  less  than 
that  of  lower  frequencies.  The  reason  for  the  smaller  delay  is  the  current  crowd- 
ing at  the  edges  of  the  intrinsic  base  region  near  the  polysilicon  contacts  which 
reduces  the  effective  current  paths  and  decreases  the  RC  delay  near  the  metal 
contact.  However,  the  delay  away  from  the  metal  contact  is  still  larger  than 
that  at  lower  frequencies  due  to  the  delay  in  FPOLY  itself. 

Though  small,  the  slight  decay  of  the  voltage  in  FPOLY  at  low  frequen- 
cies causes  asymmetric  biases  and  moves  the  minimum  voltage  point  I  toward 
FPOLY.  The  point  I  acts  as  a  current  divider  in  the  intrinsic  base  region  which 
forces  part  of  the  current  originally  flowing  through  FPOLY  at  symmetric  bias 
to  flow  through  NPOLY  Consequently,  the  current  INpoLY  from  NPOLY 
increases  with  the  decrease  of  current  IppoLY  through  FPOLY  The  contribu- 
tion of  current  from  each  base  polysilicon  contact  can  be  characterized  by  the 
effective  length  Leff  defined  in  section  4.3.3.  If  Lejf  is  longer  than  the  actual 
polysilicon  contact  length  Lactuah  the  current  contributed  from  that  contact  is 
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Figure  6.1  The  modeled  small-signal  voltage  in  FPOLY  of  a  5(im  SMB 
transistor  when  a  small-signal  voltage  with  a  magnitude  of  3mV 
is  applied  to  the  base  voltage  source.  (a)Magnitude.  (b)Phase. 
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more  than  the  current  contributed  from  the  same  contact  while  the  biases  are 
symmetric,  and  vice  versa. 

The  comparison  of  Left from  the  model  and  3D  simulations  are  plotted  in 
Fig.  6.2  as  a  function  of  frequency  and  R^.  At  low  frequencies,  Le#-of  NPOLY 
is  longer  and  Lete  of  FPOLY  is  shorter  than  Lactuai  (=5|im).  The  use  of  higher 
Rsh  increases  the  difference.  However,  the  sum  of  the  two  Le«s  approximately 
equals  two  times  Lactuai  due  to  the  division  of  the  constant  base  current 
between  the  two  polysilicon  contacts  when  the  debiasing  in  the  intrinsic  base 
region  is  negligible.  At  intermediate  frequencies,  the  magnitude  of  the  small- 
signal  voltage  in  FPOLY  decreases  slightly  while  the  phase  delay  results  in  a 
shift  from  imaginary  to  real  current.  Thus,  the  real  LeffS  become  longer  and  the 
imaginary  Le«s  are  shorter  for  both  sides  of  the  polysilicon  contacts.  As  the 
operating  frequency  further  increases,  both  the  real  and  imaginary  LeffS 
decrease  for  both  contacts  due  to  the  large  debiasing  in  FPOLY.  The  use  of 
smaller  R^  in  polysilicon  contacts  reduces  the  influences  of  debiasing  effects  at 
all  frequencies.  Therefore,  the  Leas  of  both  NPOLY  and  FPOLY  are  close  to 
^actual  compared  with  devices  having  larger  R^.  Good  matches  are  obtained 
between  the  modeled  and  simulated  Leas  at  all  frequencies  for  different  R^s 
verifying  the  accuracy  and  validity  of  the  3D  distributed  model. 

Quasi-3D  simulations  are  performed  using  the  coupling  of  2D  simula- 
tions and  Lef&  obtained  in  Fig.  6.2,  and  the  resulting  Y-qs  are  compared  with 
that  of  3D  and  2D  simulations  as  shown  in  Fig.  6.3.  For  the  DMB  structure,  Yn 
from  the  2D  simulation  is  the  same  as  Y^  of  the  3D  simulation  since  the 
voltage  distribution  in  the  polysilicon  contacts  of  the  device  are  uniform.  Yn  of 
the  SMB  structure  is  similar  to  that  of  the  DMB  at  low  frequencies.  The  use  of 
2D  simulation  can  accurately  predict  the  input  characteristics  of  the  SMB 
devices  in  this  frequency  range.  However,  Yjj  of  the  SMB  structures  deviates 
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Figure  6.3  A  comparison  of  simulated  and  modeled  Y]j  with  DMB  and  SMB 
layout  structures. 
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from  that  of  the  DMB  at  high  frequencies  due  to  the  debiasing  and  phase  delay 
in  FPOLY.  The  higher  the  sheet  resistance,  the  larger  the  debiasing  and  phase 
delay,  resulting  in  a  bigger  discrepancy  between  Yn  of  the  SMB  and  DMB 
structures.  2D  simulation  can  no  longer  properly  describe  the  input  behavior  of 
a  SMB  transistor  at  this  condition.  In  contrast,  quasi-3D  simulations  match 
the  results  of  3D  simulations  at  all  frequency  ranges  of  interest  for  various 
values  of-Rg/j.  Therefore,  the  use  of  quasi-3D  simulations  improves  the  accuracy 
of  2D  simulations  and  is  more  efficient  than  3D  simulations. 

6.3  Structure  and  Geometry  Dependent  AC  Performance 

The  distributed  effects  are  also  important  to  other  y-parameters  and 
show  different  degrees  of  influences  on  AC  behaviors  with  different  Rs^,  layout 
structures,  and  geometry  sizes.  Figure  6.4  illustrates  the  deviation  of  Y21  of  the 
SMB  devices  relative  to  Y2i  of  DMB  structure.  The  use  of  higher  Rsh  decreases 
real  Y2i  in  SMB  due  to  larger  debiasing  in  FPOLY.  Higher  Rs^  also  causes 
larger  small-signal  voltage  phase  delay  in  FPOLY  and  increases  imaginary  Y2i 
at  intermediate  frequencies.  The  imaginary  Y21  of  higher  Rsh  devices  decreases 
faster  at  higher  frequencies  due  to  the  larger  voltage  drop  in  FPOLY  following 
the  increase  in  frequency. 

The  backward  transmission  current  mainly  passes  through  the  base- 
collector  junction  capacitance  Cjbc  and  then  divides  into  two  components,  one 
flows  through  the  emitter  and  another  through  the  base.  Varying  Rs^s  in  SMB 
devices  does  not  affect  Y12  at  low  frequencies  since  only  a  small  amount  of 
current  passes  through  Cjbc.  As  the  frequency  increases,  more  current  flows 
through  Cjbc.  The  device  with  higher  Rsh  has  a  larger  effective  base  resistance 
which  reduces  the  current  passing  through  it.  Consequently,  Y12  of  those 
devices  becomes  smaller  at  high  frequencies  as  shown  in  Fig.  6.5.  Larger  effec- 
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Figure  6.4  The  Y21  plot  of  DMB  and  SMB  structures  with  different  Rsh 
obtained  from  3D  simulations. 
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tive  base  resistance  also  increases  the  potential  in  the  base  region  forcing  more 
electrons  to  be  injected  from  the  emitter  into  the  collector.  Thus,  Y22s  of  the 
devices  with  larger  Rsh  are  larger  than  that  of  smaller  ones  at  low  and  inter- 
mediate frequencies  as  illustrated  in  Fig.  6.6.  The  collector  current  becomes 
dominant  by  the  current  flowing  through  Cjbc  as  the  frequency  increases.  The 
larger  effective  base  resistance  in  SMB  structures  induce  less  phase  angle 
which  increases  the  real  and  decreases  the  imaginary  collector  currents,  and 
likewise  the  real  and  imaginary  part  of  Y22.  Though  the  use  of  different  values 
o£Rsh  significantly  changes  the  y-parameters,  the  cut-off  frequencies  fy  of  these 
devices  (~32GHz)  are  the  same  since  fT  is  mainly  determined  by  the  properties 
of  the  base-emitter  junction  and  the  intrinsic  base  regions. 

The  distributed  effects  acting  in  FPOLY  of  the  SMB  structure  strongly 
influence  the  small-signal  current  behavior  at  high  frequencies.  For  the  inter- 
digitated  structure,  similar  distributed  effects  take  place  in  both  base 
polysilicon  contacts.  The  modulation  of  current  at  high  frequencies  is  more 
significant  than  in  the  SMB  devices  since  debiasing  occurs  in  both  fingers  of 
interdigitated  structure.  Figure  6.7  compares  Yn  for  the  DMB,  SMB  and  inter- 
digitated  structures  with  i?s^=12.5Q/sq.  It  is  seen  that  both  the  real  and 
imaginary  parts  of  Y^  of  the  interdigitated  structure  are  smaller  than  that  in 
the  SMB  device  at  high  frequencies,  and  the  deviation  from  Yn  of  DMB  is  more 
significant  compared  with  SMB  devices.  Thus,  the  characterization  of  the  inter- 
digitated structure  with  2D  simulations  results  in  a  more  significant  error  than 
the  SMB  structure. 

The  input  and  output  small-signal  current  driving  capability  is  usually 
assumed  to  be  linearly  proportional  to  the  stripe  length,  and  the  phase  delay  in 
the  polysilicon  contacts  are  neglected.  For  the  DMB  structure,  there  are  no 
distributed  effects  in  the  polysilicon  contacts.  Therefore,  |Yn|  and  |Y2i|  are 
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Figure  6.5  The  Y12  plot  of  DMB  and  SMB  structures  with  different  Rsh 
obtained  from  3D  simulations. 
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Figure  6.7  A  comparison  of  Yn  between  DMB,   SMB   and  interdigitated 
structures  withi?s^=12.5Q/sq  obtained  from  3D  simulations. 
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linearly  proportional  to  the  stripe  length  as  shown  in  Fig.  6. 8. (a)  and  (b), 
respectively.  The  phase  delays  are  also  the  same  regardless  of  the  length  of  the 
stripe  used  as  illustrated  in  Fig.  6.9.  However,  in  the  presence  of  distributed 
effects,  debiasing  along  FPOLY  in  the  SMB  and  the  fingers  in  interdigitated 
structures  reduces  the  input  and  output  current  drive  and  degrades  the  linear- 
ity at  high  frequencies.  The  signal  propagating  along  FPOLY  and  the  fingers 
also  requires  a  longer  time  compared  with  that  of  the  DMB  structure  resulting 
in  additional  phase  delay.  A  longer  stripe  further  degrades  the  current  drive 
capability  and  enlarges  the  RC  delay  unless  a  very  low  value  ofRsh  is  utilized. 
The  sheet  resistance  of  the  base  polysilicon  contact  plays  a  key  role  in  deter- 
mining current  driving  capability  and  phase  delay  in  SMB  and  interdigitated 
devices.  For  example,  the  input  current  driving  capability  of  a  lOfim  SMB 
device  with  i?s/j=12.5Q/sq  is  18%  less  and  phase  delay  is  12°  more  than  a  DMB 
structure.  If  ai?s^  of  3.1£2/sq  is  used,  the  current  loss  and  additional  RC  delay 
reduces  to  9.3%  and  5°,  respectively.  Thus,  a  good  silicidation  in  the  base 
polysilicon  contacts  can  significantly  improve  the  characteristics  of  SMB  and 
interdigitated  BJTs. 

The  interdigitated  structure  has  worse  current  driving  capability  and 
larger  phase  delay  than  the  SMB  device  since  the  debiasing  effects  and  RC 
delay  take  place  in  both  polysilicon  contacts.  Due  to  distributed  effects  in  the 
base  polysilicon,  at  its  worst,  the  current  in  a  SMB  device  is  totally  supplied  by 
NPOLY.  For  the  interdigitated  structure,  the  distributed  effects  keep  affecting 
the  current  and  phase  delay  as  the  frequencies  and  stripe  length  increase. 
Consequently,  the  input  and  output  current  drive  of  SMB  devices  increase 
monotonically,  while  the  current  supply  of  interdigitated  devices  tends  to 
saturate  as  the  stripe  length  increases.  At  the  same  time,  the  phase  delay  of 
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the  SMB  is  limited  while  that  of  the  interdigitated  structure  keeps  increasing 
following  the  increase  of  stripe  length. 

6.4.  The  Importance  of  Distributed  Effects  on  Circuit  Design 

The  use  of  different  layout  structures  does  not  affect  the  small-signal 
current  gain  and  cut-off  frequency.  However,  there  is  a  significant  difference  of 
Yji  and  a  minor  difference  in  other  y-parameters  between  various  layout  struc- 
tures. The  discrepancy  requires  accordant  design  of  matching  networks  in  an 
amplifier  to  achieve  stability  and  power  gain  criterion. 

Figure  6.10  illustrates  the  definition  of  reflection  coefficients  of  a  two- 
port  network.  In  terms  of  reflection  coefficients,  the  conditions  for  uncondi- 
tional stability  of  the  network  at  a  given  frequency  are  [Gon84] 

|rs|<i,  (6-i) 

|rL|  <  i ,  (6-2) 

<  1 ,  (6-3) 


Fin 


Q       ,  S\2S2\TL 


and 


TOOT1 


1  ~S22TL 


Q       ,  ^12^21  rS 


<  1 .  (6-4) 


'Ul  s 

There  are  many  sets  of  Ts  and  Fl  that  can  produce  stable  operation.  Figure 
6.11  plots  the  stability  circle  of  T-^  in  the  Smith  chart  for  different  layout  struc- 
tures. The  area  inside  the  circles  indicates  the  value  of  Tl  producing  unstable 
operation.  It  is  found  that  this  device  is  unconditionally  stable  for  the  SMB  and 
interdigitated  structures  with  any  passive  load  and  tends  to  be  unconditionally 
stable  for  the  DMB  device  unless  real  TL  approaches  one.  Any  passive  input 


146 


VsZO°Q 


rs      riN 


rouT    rL 


Zs 


Two-port 
network 


-IN 


-OUT 


Figure  6.10  The  definition  of  reflection  coefficients  of  a  two-port  network.  (After 
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Figure  6.11  rL  stability  circles  and  constant  power  gain  circles  in  the  Smith 
chart  with  three  layout  structures.  The  specified  power  gain  Gp  is 
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impedance  in  the  network  also  produces  an  unconditional  stability  for  the  SMB 
and  interdigitated  structures  and  tends  to  be  unconditionally  stable  for  the 
DMB  device  unless  |  Eg  |  approaches  one.  Although  all  the  structures  tend  to 
be  unconditionally  stable  in  this  case,  it  is  seen  that  the  use  of  various  layouts 
corresponds  to  different  stability  circles.  Thus,  the  available  Tg  and  T^  can  be 
different  for  different  layout  structures  in  some  other  devices. 

In  microwave  circuit  design,  power  gain  is  one  of  the  most  important 
figures  of  merit,  and  the  operating  power  gain  Gp  is  independent  of  the  source 
impedance.  Therefore,  a  design  procedure  based  on  Gp  for  both  unconditionally 
stable  and  potentially  unstable  transistors  is  simple  and  commonly  used  for 
practical  design  [Gon84] .  Gp  is  defined  as  the  ratio  of  power  delivered  to  the 
load  vs.  power  input  to  the  network  and  is  given  by 
(  \ 
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A  -  SuS22-SnS2l. 


(6-6) 


Both  input  and  output  characteristics  of  the  transistor  affect  Gp.  Therefore,  Gp 
also  depends  on  the  layout  structures  used.  A  power  gain  circle  in  a  T^  plane  of 
the  Smith  chart  gives  the  available  region  of  T^  for  which  a  specific  gain  can  be 
achieved  at  a  given  frequency.  The  power  gain  circles  with  Gp-40  at  a 
frequency  of  10GHz  for  different  layout  structures  are  plotted  in  Fig.  6.11 
together  with  the  output  stability  circles.  For  both  the  5  and  10  |im  devices,  the 
DMB  structure  allows  the  widest  range  of  T^  that  can  achieve  the  specified  Gp. 
There  are  more  constraints  for  available  T^  of  the  SMB  device  compared  with 
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that  of  the  DMB,  and  the  interdigitated  structure  has  the  smallest  range  of 
available  T^. 

As  mentioned  in  Chapter  3  and  4,  2D  device  simulations  can  only 
describe  devices  similar  to  the  DMB  structure  with  uniform  distributed  voltage 
and  current  in  the  z-direction.  However,  2D  simulation  is  still  one  of  the  main 
tools  in  device  design  independent  of  the  layout  structures.  This  means  the  3D 
distributed  effects  in  the  device  and  circuit  design  are  neglected.  So,  how  much 
difference  will  it  make  if  these  effects  are  not  taken  into  account?  Assuming  an 
output  matching  network  is  designed  according  to  2D  simulations  (DMB  struc- 
ture) and  the  goal  is  to  achieve  a  power  gain  of  40,  any  T-^  in  the  gain  circle  of 
DMB  can  be  used  except  for  real  TL  ~  1  because  of  stability  considerations.  If  a 
TL  is  chosen  from  point  A  indicated  in  Fig.  6.10,  the  output  matching  network 
is  simply  a  resistor  of  433Q  and  an  inductor  of  1.98nH  in  series  for  a  5|im  stripe 
BJT.  If  this  matching  network  is  applied  to  a  SMB  device,  a  1.21dB  power  gain 
loss  results  and  a  1.75dB  loss  occurs  in  an  interdigitated  structure.  The  degra- 
dation is  even  worse  for  longer  stripe  transistors.  Applying  the  matching 
network  designed  for  the  DMB  device  to  a  SMB  device  causes  a  2.39dB  power 
gain  loss  and  a  loss  of  4.44dB  in  the  interdigitated  transistor.  The  power  gain 
for  the  BJT  with  different  layouts  are  summarized  in  Table  6.1  and  6.2. 

6.5.  Summary 

The  accuracy  of  the  AC  distributed  model  derived  in  Chapter  4  and  the 
quasi-3D  simulations  have  been  verified  by  3D  device  simulations  in  this 
chapter.  The  device  characteristics  dependencies  on  the  layout  structure  and 
sheet  resistance  due  to  the  presence  of  3D  distributed  effects  are  also  studies. 
It  is  shown  that  the  use  of  different  types  of  layout  structures  and  geometries 
result  in  a  significant  difference  in  the  input  characteristics  of  the  same 
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transistor  unless  a  very  low  value  of  sheet  resistance  is  used  in  the  base  polysil- 
icon.  A  good  silicided  polysilicon  contact  greatly  eliminates  the  debiasing  and 
improves  the  small-signal  current  driving  capability  for  long  stripe  BJTs. 

The  varied  characteristics  due  to  the  use  of  different  layout  structures 
require  individual  design  of  the  matching  network.  It  is  proved  that  applying 
an  output  matching  network  designed  for  a  DMB  structure  to  a  SMB  and  inter- 
digitated  structure  results  in  a  significant  degradation  in  the  power  gain.  The 
stability  properties  can  also  be  degraded  in  other  cases.  Therefore,  the  distrib- 
uted effects  must  be  taken  into  account  in  the  design  of  SMB  and  interdigitated 
devices. 

2D  device  simulations  fail  to  account  for  the  voltage  and  current  varia- 
tion in  the  base  polysilicon  region  resulting  in  a  significant  error  in 
determining  AC  characteristics  for  SMB  and  interdigitated  structures,  while 
3D  device  simulations  demand  very  extensive  computer  resources.  An  AC 
FLOODS  simulation  with  2100  nodes  requires  1.5G  of  memory  space  running 
on  a  Sun  Ultrasparc  2  computer  for  more  than  30  hours.  In  contrast,  the  quasi- 
3D  simulations  predict  the  AC  characteristics  accurately  and  demand  much 
less  computer  memory  and  simulation  time  making  it  is  useful  for  practical 
device  and  circuit  design. 
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Table  6.1.  The  power  gain  of  a  5|Jjn  transistor  with  different  layouts 


Layout 


DMB 


SMB 


Interdigitated 


Power  Gain 


40.0  (16.02dB) 


30.0  (14.77dB) 


26.7  (14.27dB) 


Table  6.2.  The  power  gain  of  a  10(im  transistor  with  different  layouts 


Layout 


DMB 


SMB 


Interdigitated 


Power  Gain 


40.0  (16.02dB) 


23.0(13.62dB) 


14.4(11.58dB) 


CHAPTER  7 

THE  DISTRIBUTED  EFFECTS  ON  THE  AC  AND  TRANSIENT 

PERFORMANCE  OF  A  0.18um  MOSFET 


7.1  Introduction 

MOSFETs  are  the  dominant  devices  in  current  digit  circuit  applications, 
and  the  applications  are  also  being  extended  to  radio  frequency  applications  as 
the  feature  sizes  are  reduced.  The  conventional  layout  and  cross  section  of  a 
MOSFET  with  a  polysilicon  gate  is  illustrated  in  Fig.  7.1.  The  metal  contacts 
of  a  gate  are  located  at  the  ends  of  the  polysilicon.  The  input  signal  needs  to 
propagate  from  metal  contacts  through  the  gate  to  control  the  charges  in  the 
channel.  At  high  frequency  operation,  an  input  signal  RC  delay  takes  place 
along  the  gate  due  to  the  distributed  gate  resistance  and  capacitance.  The 
scaled  geometries  in  the  polysilicon  gate  length  (channel  length),  thickness, 
and  the  gate  oxide  thickness  increase  the  resistance  and  capacitance  resulting 
in  a  larger  RC  delay.  This  delay  becomes  important  if  the  transistor  width  is 
increased  to  obtain  high  current  drive  and  large  transconductance.  This  effect 
is  especially  noticeable  in  line  drivers  and  output  buffers  used  in  digital  circuits 
and  low  noise,  high  gain  amplifiers  employed  in  analog  applications  [Raz94]  in 
which  a  wide  gate  is  usually  required. 

The  delay  time  is  one  of  the  figures  of  merit  in  measuring  the  perfor- 
mance of  a  MOSFET.  A  low  gate  resistance  using  silicided  technologies  is 
usually  emphasized  in  device  design  to  avoid  the  RC  delay  due  to  the  distrib- 
uted gate  resistance  [Mii94,  Raz94].  However,  the  relation  between  the  RC 
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Figure  7.1.  The  layout  and  cross  section  of  a  conventional  nMOSFET. 
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delay,  device  performance,  and  the  gate  sheet  resistance  at  different  frequen- 
cies is  not  clarified. 

In  this  chapter,  the  implementation  of  the  transient  simulation  scheme 
in  FLOODS  is  briefly  discussed.  The  RC  delay  in  the  transient  and  AC  opera- 
tions of  a  0.18(im  MOSFET  with  various  gate  sheet  resistances  is  examined 
based  on  3D  FLOODS  simulations. 

7.2  Transient  Simulations 
The  general  form  of  the  semiconductor  equations  is  given  by 


~q(z(t))  +  f(t,z(t))  =  0 


(7-1) 


where  z(t)  is  the  solution  of  the  equations  in  nodes  at  time  t ,  q(z)  is  the  transient 
and  f(t,z)  is  the  steady-state  term  in  the  equations.  The  z(t)s  and  q(z)s  for  each 
semiconductor  equation  used  in  FLOODS  are  summarized  in  Table  7.1  where 
Ec  and  Ev  are  the  conduction  and  valance  bands,  respectively,  kg  is  the  Boltz- 
mann  constant,  C^  is  the  lattice  specific  heat,  and  y  is  the  Fermi  degenerate 
coefficient  (refer  to  p.  24,  [Lia94c]).  The  TR-BDF2  method  proposed  by  Bank 
[Ban85]  is  used  in  FLOODS  to  discretize  Eq.  (7-1)  in  time  domain. 

Table  7.1 


Equation 

z(t) 

q(z) 

Poisson 

potential  \\f 

0 

Electron  Continuity 

electron  n 

n 

Hole  Continuity 

hole/? 

P 

Electron  Energy  Balance 

electron  temperature  Tn 

n(Ec+1.5y/c57>i) 

Hole  Energy  Balance 

hole  temperature  T 

p(-Ev+l.  SykgTp) 

Lattice  Self-Heating 

lattice  temperature  TL 

CJl 
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The  total  current  Jtot  in  the  terminals  of  semiconductor  devices  is  given 
by 

Jtot   =   Jn+Jp+Jdlsp  (7-2) 

where  Jn  and  Jp  are  the  electron  and  hole  conduction  currents,  respectively, 
and  Jdisp  =  e-j-  is  the  displacement  current  due  to  the  presence  of  space 
charges  at  the  terminals,  where  E  is  the  electric  field  and  E  is  the  permittivity 
in  the  semiconductor  material.  The  space  charges  in  the  terminals  disappear 
fast  due  to  the  very  short  dielectric  relaxation  time.  Therefore,  the  displace- 
ment current  is  usually  not  important  in  transient  operations  of  contemporary 
semiconductor  devices.  Figure  7.2  shows  the  transient  response  of  a  diode 
obtained  from  FLOODS  simulations.  It  is  seen  that  the  displacement  current 
Jdisp  is  only  important  for  the  first  30  femto  seconds,  then  the  carrier  conduc- 
tion currents  Jn  and  Jp  dominate  the  total  current  for  the  rest  of  operations. 
Although  J  (lisp  is  n°t  important  in  the  time  domain  of  interest,  it  is  found  that 
the  carriers  in  the  semiconductor  are  forced  to  response  abnormally  fast  iiJ^isp 
is  absent  from  the  total  current  in  Eq.  (7-2),  which  results  in  numerical  insta- 
bility and  significant  error  in  the  solutions  of  the  device  simulations. 

7.3  Distributed  Effects  on  the  Transient  Performance  of  MOSFETs 
The  simulated  structure  of  a  0.18|im  nMOSFET  with  a  gate  oxide  thick- 
ness of  40  A  is  illustrated  in  Fig.  7.3.  The  width  of  this  transistor  is  10|im  and 
there  are  2316  nodes  in  the  simulation  structure.  To  properly  represent  the 
location  of  metal  contact  in  a  real  device  structure,  the  gate  metal  contact  in 
this  simulation  structure  is  attached  to  the  face  A.  The  input  signal  needs  to 
propagate  from  the  metal  contact  through  the  polysilicon  gate  to  control  the 
output  characteristics  in  this  Edge  Contact  (EC)  structure  and  a  RC  delay 
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Figure  7.2.  The  transient  current  of  a  diode  when  a  ramp  voltage  is  applied  at 
£=0.  The  inserted  picture  illustrates  the  circuit  diagram  of  the 
simulations. 
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Figure  7.3.  The  simulation  structure  of  a  0.18|im  MOSFET.  The  material  of  the 
bottom  cube  is  silicon  and  the  cube  on  the  top  is  polysilicon  gate. 
There  is  a  40  A    gate  oxide  between  silicon  and  gate  polysilicon. 
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Figure  7.4.  The  transient  source  voltage  in  the  turning  on  of  a  DRAM  obtained 
from  3D  FLOODS  simulations.  The  inserted  picture  illustrates  the 
circuit  diagram  of  the  simulations. 
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takes  place  along  the  signal  path.  The  sheet  resistance  of  the  polysilicon  gate 
is  controlled  by  using  different  mobility  values  in  the  gate  material  which  is 
independent  of  the  mobilities  in  other  materials.  For  comparison  reasons,  an 
artificial  device  in  which  the  gate  metal  contact  is  attached  to  the  Face  B  is  also 
used  in  the  study.  Since  the  metal  contact  covers  the  entire  polysilicon  gate  top 
surface,  there  is  no  RC  delay  along  the  gate  in  this  Full  Contact  (FC)  structure. 

The  dynamic  random  access  memory  (DRAM)  cell  contains  a  capacitor 
to  store  the  charge  and  a  MOSFET  switch  to  read  and  write  the  stored  charge. 
The  performance  of  a  DRAM  cell  depends  on  the  switching  speed  of  the 
MOSFET  and  the  time  required  for  the  MOSFET  to  charge/discharge  the 
capacitance.  A  delay  in  the  input  signal  slows  down  the  turning  on  and  off  of 
the  transistor  and  increases  the  charge/discharge  time.  The  transient  charac- 
teristics of  the  turning  on  of  a  DRAM  cell  with  the  0.18um  MOSFET  are  plotted 
in  Fig.  7.4.  In  the  simulations,  the  gate  sheet  resistance  of  the  FC  structure  is 
2Q/sq  and  that  of  the  EC  structure  is  570/sq.  A  step  voltage  of  2V  is  applied  to 
the  gate  metal  contact  at  t=0.  The  response  of  the  source  terminal  voltage  is 
seen  to  significantly  lag  behind  to  the  input  voltage  for  both  the  FC  and  EC 
structures.  However,  the  difference  of  the  voltage  response  between  FC  (no  gate 
delay)  and  EC  (with  gate  delay)  devices  is  small  implying  the  RC  delay  in  the 
gate  is  insignificant  compared  with  the  DRAM  capacitance  charging  time. 

The  gate  sheet  resistance  is  usually  lower  than  57H/sq  for  contemporary 
MOSFET  technologies,  and  the  transistor  width  is  normally  not  this  wide  for 
most  of  digital  applications.  Therefore,  the  distributed  effects  can  usually  be 
neglected  in  the  MOSFETs  for  digital  applications. 
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7.4  Distributed  Effects  on  the  AC  Performance  of  MOSFETs 
The  applications  of  MOSFETs  are  being  extending  to  the  radio 
frequency  applications  as  the  feature  sizes  are  reduced.  A  cut-off  frequency  of 
93GHz  has  been  reported  for  a  0.08um  channel  MOSFET  [Mii94].  The  gate 
small-signal  current  increases  as  the  operating  frequency  increases  due  to  the 
larger  displacement  current  passing  through  the  gate  oxide.  The  transistor 
width  in  the  high  frequency  applications  is  usually  large  to  achieve  high 
current  drive  and  large  transconductance.  Therefore,  the  AC  debiasing  and 
phase  delay  due  to  the  distributed  gate  resistance  and  capacitance  become 
important  at  high  frequency  applications. 

|  Yn  |  of  a  MOSFET  increases  following  the  increasing  frequency  due  to 
the  larger  small-signal  current  flowing  through  the  gate  at  higher  frequencies. 
The  phase  delay  of  Yn  is  also  larger  because  of  the  more  significant  RC  delay 
along  the  gate  current  path.  The  use  of  a  larger  gate  sheet  resistance  increases 
the  gate  impedance  resulting  in  a  smaller  |  Yn  |  and  a  more  significant  phase 
delay  in  Yn.  The  simulated  Yn  of  the  MOSFET  is  shown  in  Fig.  7.5  as  a 
function  of  frequency  and  gate  sheet  resistance  at  Vqs=3V  and  Vrjs=2V.  There 
is  little  phase  delay  for  FC  structure  at  high  frequencies  since  there  is  no  RC 
delay  in  the  gate.  Compared  to  the  FC  structure,  |Yn|  of  EC  structure  is 
smaller  since  the  effective  gate  impedance  increases  due  to  the  longer  current 
path.  The  phase  delay  of  Yn  in  EC  device  is  also  lagged.  The  difference  of  Yn 
between  the  FC  and  EC  structures  with  a  Rsh=57Q/sq  is  significant  at  high 
frequencies.  If  Rsh  is  reduced,  the  difference  becomes  smaller  and  in  fact  negli- 
gible if  Rsh  is  less  than  lOQ/sq. 

As  the  operating  frequency  increases,  the  increasing  small-signal 
current  causes  a  larger  debiasing  along  the  gate  which  reduces  the  voltage 
across  the  gate  and  decreases  the  drain  current.  The  larger  input  signal  delay 
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in  the  gate  also  lags  the  response  of  the  output  drain  current  at  higher  frequen- 
cies. Thus,  |  Y2i  |  decreases  and  phase  delay  of  Y2i  increases  following  the 
increase  of  operating  frequency.  Figure  7.6  shows  the  transconductance  Y2i  at 
various  frequencies  and  Rshs.  |  Y2i  |  of  FC  device  is  independent  of  the 
frequency  since  there  is  no  debiasing  along  the  gate.  A  reduction  of  |  Y21  |  is 
observed  at  high  frequencies  for  the  EC  structures.  Larger  Rsh  causes  more 
significant  decrease  due  to  a  bigger  small-signal  ohmic  voltage  drop  along  the 
gate.  The  phase  delay  of  Y21  for  the  EC  devices  is  also  larger  at  high  frequen- 
cies compared  with  that  of  the  FC  device  due  to  the  additional  gate  RC  delay. 
However,  the  difference  between  the  FC  and  EC  structures  is  insignificant  if 
the  gate  sheet  resistance  is  smaller  than  lOQ/sq.  The  analysis  on  the  Y12  and 
Y22  also  reveals  the  distributed  effect  degrades  the  characteristics  of  devices  at 
high  frequencies  for  large  gate  sheet  resistance  and  becomes  negligible  if  Rsh 
is  less  than  lOQ/sq. 

Therefore,  a  MOSFET  with  a  large  gate  sheet  resistance  suffers  from  the 
distributed  effect  at  high  frequencies  and  the  performance  degrades.  The 
silicided  gate  technologies  eliminate  the  distributed  effect  and  improve  the 
performance.  If  the  gate  sheet  resistance  is  less  than  lOQ/sq,  this  effect  can  be 
neglected  in  the  device  and  circuit  design. 

7.4  Summary 
The  distributed  effects  on  the  performance  of  an  advanced  0.18|im 
nMOSFET  have  been  studied  based  on  3D  device  simulations.  It  has  been 
shown  that  the  RC  delay  due  to  the  distributed  gate  resistance  and  capacitance 
contributes  little  to  the  delay  time  of  a  DRAM  cell  even  if  a  large  gate  sheet 
resistance  is  used.  Therefore,  the  performance  of  MOSFETs  in  digital  applica- 
tions is  fairly  insensitive  to  silicided  gate  technologies.  The  distributed  effects 
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are  significant  for  MOSFETs  with  large  gate  sheet  resistance  in  radio 
frequency  applications.  The  use  of  silicided  gate  technologies  is  necessary  to 
improve  the  performance  of  high  speed  MOSFETs.  However,  if  a  reasonably 
low  polysilicon  gate  sheet  resistance  is  utilized,  the  distributed  effects  and  RC 
delay  in  the  gate  can  be  neglected  in  the  radio  frequency  device  and  circuit 
design. 
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Figure  7.5.  The  simulated  Yn  as  a  function  of  frequency  and  Rsh. 
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Figure  7.6.  The  simulated  Y2i  as  a  function  of  frequency  and  Rs^. 


CHAPTER  8 
SUMMARY  AND  FUTURE  WORK 


8.1  Summary 

The  voltage  and  current  in  the  base  region  of  a  BJT  and  in  the  polysili- 
con  gate  region  of  a  MOSFET  exhibit  the  distributed  effects.  This  distributed 
effect  is  significant  when  a  device  operates  at  high  levels  of  current  injection  or 
at  high  speed.  Doing  device  and  circuit  design  without  taking  into  account  this 
effect  can  result  in  degraded  device  and  circuit  performance. 

Distributed  effects  in  the  polysilicon  and  intrinsic  base  regions  of  stripe 
B JTs  with  double  polysilicon  technology  were  studied  in  detail  in  this  disserta- 
tion. Analytical  models  describing  the  DC  and  AC  voltage  and  current 
distributions  in  the  devices  were  presented  in  Chapters  3  and  4,  respectively. 
This  quasi-3D  model  properly  predicted  the  input  characteristics  of  a  single 
metal  base  contact  BJT  with  double  polysilicon  technologies.  A  quasi-3D  device 
simulation  method  was  also  demonstrated  in  chapter  4  to  accurately  describe 
the  3D  effects  using  2D  device  simulations,  which  enables  engineers  to  design 
devices  and  circuits  efficiently. 

The  influences  of  the  distributed  effect  on  various  types  of  BJT  layout 
structures,  geometry  sizes,  and  polysilicon  contact  resistances  were  studied  in 
Chapters  5  and  6  for  DC  and  AC  operations,  respectively.  It  was  shown  that 
very  low  resistance  polysilicon  is  necessary  for  high  performance  BJTs.  There- 
fore, a  silicided  polysilicon  process  is  one  of  the  key  factors  in  determining  the 
device  performance  of  a  BJT.  The  use  of  different  layout  structures  and  geome- 
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try  sizes  also  results  in  significant  differences  in  input  and  output 
characteristics  of  BJTs.  The  linear  relation  between  current  drive  and  geome- 
try sizes  is  invalid  in  devices  with  long  stripe  operating  at  high  biases  or  high 
frequencies.  A  design  without  taking  into  account  the  distributed  effect  was 
shown  to  result  in  a  poor  circuit  performance. 

The  distributed  effect  was  also  observed  in  an  advanced  MOSFET  This 
effect  significantly  degrades  the  transconductance  and  delay  the  output 
response  if  a  large  gate  sheet  resistance  is  used  for  high  frequency  applications. 
The  silicided  polysilicon  technologies  are  necessary  to  improve  the  performance 
of  high  speed  MOSFETs.  However,  the  effect  is  less  important  on  the  high 
frequency  performance  of  MOSFETs  if  an  adequate  polysilicon  gate  sheet 
resistance  is  utilized.  The  results  also  showed  the  distributed  effect  does  not 
significantly  affect  the  MOSFET  performance  in  digital  applications  even  if  a 
large  gate  sheet  resistance  is  utilized. 

8.2  Future  Work 
The  following  subjects  are  suggested  to  continue  the  applications  of 
distributed  effects: 

(1)  Implementing  the  automatic  quasi-3D  simulation  capability  in  FLOODS. 
Once  the  3D  device  information  are  given,  all  the  device  parameters  and  coeffi- 
cients in  the  model  can  be  determined  from  2D  device  simulations.  The 
effective  lengths  are  then  computed  from  the  analytical  equations  and  feed  to 
the  simulator  for  quasi-3D  simulations.  This  procedure  can  be  included  in  a 
macro  function  of  FLOODS  to  automatically  perform  the  quasi-3D  simulations. 

(2)  Including  the  automatic  extraction  function  for  base  resistance  in  FLOODS. 
Base  resistance  is  one  of  the  most  important  and  most  difficult  parameters  in 
BJTs  to  extract.  All  the  device  parameters  and  coefficients  required  in  the 
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model  to  compute  the  base  resistance  can  be  obtained  from  2D  device  simula- 
tions. Therefore,  a  post-processor  reads  the  2D  simulation  results  can  compute 
the  base  resistance  directly. 

(3)  Including  the  distributed  effects  in  a  compact  model.  A  circuit  simulator  is 
the  main  tool  in  current  integrated  circuit  design.  The  accuracy  of  the  simula- 
tion results  depends  on  the  accuracy  of  the  compact  device  models  used  in  the 
simulator.  AC  base  resistance  is  one  of  the  most  important  parameters  in  the 
compact  model.  The  impedance  semicircle  method  based  on  the  simple  lumped 
circuit  model  [Fus95]  is  one  way  to  extract  the  AC  base  resistance.  Figure  8.1 
shows  an  impedance  semicircle  in  a  complex  plane  obtained  from  3D  simula- 
tions with  a  5|_im  stripe  length  SMB  BJT  similar  to  the  structure  used  in 
chapter  5.  The  base  resistance  is  extracted  from  the  intersection  of  the  semicir- 
cle and  the  real  axis  at  high  frequencies.  However,  the  input  impedance 
obviously  cannot  be  expressed  by  a  circle  equation  due  to  the  input  distributed 
resistance  and  capacitance.  Therefore,  the  inclusion  of  this  effect  in  the 
compact  models  is  important  for  high  speed  circuit  designs. 
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Figure  8.1.  Locus  of  input  impedance  in  complex  plane. 
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